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Abstract 



In this scries of lectures a method is developed to compute one-loop shifts to classical masses 
of kinks, multi-component kinks, and self-dual vortices. Canonical quantization is used to show 
that the mass shift induced by one-loop quantum fluctuations is the trace of the square root of the 
differential operator governing these fluctuations. Standard mathematical techniques are used to 
deal with some powers of pscudo-diffcrcntial operators. Ultraviolet divergences arc tamed by using 
generalized zeta function regularization methods and, then performing zero-point energy and mass 
renormalizations. Information about the meromorphic structure of the generalized zeta function 
of the second-order fluctuation operator K around the classical solution is obtained from the K- 
heat equation kernel via the Mellin transform. In particular, the high-temperature expansion 
of the partition function provides the residua at the poles of the generalized zeta function in 
k> , terms of the Seeley coefficients of the asymptotic approximation. In this way a formula is derived 

; I ' that allows computation of one-loop mass shifts for kinks, multi-component kinks, and self-dual 

I Abrikosov-Niclsen-Olesen vortices. Numerical results for the Seeley coefficients as well as the mass 

shifts, obtained by means of a Mathematica environment implemented on a standard PC, are 
offered. A qualitative analysis of the outcome shows a common trend in the mass shift of the 
three types of topological defects analyzed. A comparison with exact results is presented whenever 
possible, i.e., for the kink and the TA'l kink, respectively, of the Xcj)'^ and BNRT models. One- 
loop renormalization of the planar Abelian Higgs model requires use of the Feynman-'t Hooft 
renormalizablc gauge, in the vacuum sector, or the background gauge, in vortex sectors. Faddccv- 
Popov ghosts that restore unitarity are dealt with in the Hamiltonian framework in a novel fashion. 
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1 Introduction 



1.1 A brief history of soliton quantization 

We start this long Introduction by offering a sliort and biased fiistory of classical solitons and their 
quantization. The emphasis will be oriented towards the topics to be analyzed in this set of Lectures, 
skipping many important aspects of such a broad and fertile subject. 



• Solitons and solitary waves 

Traditionally, wave phenomena in nature have been distinguished by their dispersive character, 
i.e., the property by which propagating waves eventually fade away in finite time. Fascination 
with the soliton phenomenon started with the "experimental" observation of the Scottish engi- 
neer Scott-Russell circa 1870 in a Edinburgh channel:"^ solitary wave travels without changing 
its shape, size, or, speed\ [I]. 

Linear wave equations only admit traveling or solitary wave solutions if the dispersion law linking 
the frequency of the wave motion with the wave vector of a "monochromatic" component is 
linear, because in such a case all the waves in a wave packet travel with the same speed without 
interferences between them. PDEs of this type are very rare but very well known: they are 
essentially variations of the free-string and massless Dirac equations. Thus, the impact of the 
discovery by Korteweg-de Vries around 1905 of their non-linear PDE describing wave motion 
in shallow waters in channels has been enormous. Besides providing a mechanism by which 
the non-linearity balances the dispersive character of the KdV equation, circa 1965 Kruskal, 
Miura, Lax and others showed that that this magic equation can be completely solved despite 
its complexity. Between the solutions of the KdV equation there are solitary waves that keep 
their shape, height and speed during the propagation, thus providing a theoretical explanation 
for Scott-Russell traveling lumps of water. Due to complete integrability, KdV solitary wave 
solutions not only keep their shapes in free propagation but also survive collisions without 
damage, just as fundamental particles survive scattering at not too high energies. For this 
reason, the non-dispersive solutions of the KdV equation were christened as solitons. 

In fact, the solution of the KdV equation was the main impulse that led to the creation and 
development of new ideas and techniques of extraordinary importance in Mathematical Physics 
over the last fifty years, such as the inverse scattering method (Kruskal/Miura), Lax pairs and 
non-linear compatibility conditions (P. Lax), classical spectral transforms (Sakharov), etcetera. 
Also, old and almost forgotten methods such as the Backlund transformation (W. Lamb) or 
highly sophisticated ideas of algebraic geometry (Novikov/Dubrovine) found a new playground 
for application. Most remarkably, similar unexpected properties were discovered in other non- 
linear PDE, such as the non-linear Schrodinger equation and the sine-Gordon equation. Amaz- 
ingly, both equations govern the dynamics of real physical systems (thus displaying the soliton 
phenomenon): the non- linear Schrodinger equation governs some phenomena in non-linear op- 
tics; the sine-Gordon equation (discovered in geometry) explains the Josephson effect in semi- 
conductor physics, etcetera . 

• Topological defects in condensed matter physics and cosmology 

In several of these systems and in other higher-dimensional relatives, there are topological reasons 
for the strong stability of solitary waves. By this statement we mean that non- linear PDE equa- 
tions of this type are sometimes the variational Euler-Lagrange equations of some Lagrangian 
functional and, in fewer cases that are essentially one-dimensional, also admit a Hamiltonian 
formulation amenable to a sum of infinite angle-action variables. As a general feature, the con- 
figuration space is the sum of several (frequently infinite) topologically disconnected sub-spaces. 
Because temporal evolution is a homotopy transformation, field configurations in different topo- 
logical sectors cannot evolve into each other. For this reason, lumps arising as absolute minima 
of the energy in each topological sector are called topological defects. In a more complex physical 
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system, in superconductors of Type II (some alloys below the critical temperature) magnetic flux 
tubes were discovered by Abrikosov circa 1957 and were understood by him to be topological 
defects arising in the Ginzburg-Landau phenomenological theory of superconductivity. 

Similar topological defects forming tubes along a central line were also discovered in liquid 
crystals and quantum fluids by other Nobel laureates such as de Gennes and Legget, respectively 
in 1973 and 1978, who also found domain walls, point defects and textures in these exotic 
materials. The topological and group theoretical roots of these extended structures arising in 
nematic and cholesteric liquid crystals or in phases A and B of helium 3 have been studied in 
depth by Mermin, Michel and others. 

More recently, Kibble and others, circa 1989, studied how Cosmology would be affected by the 
existence of domain walls in the Universe itself. Following this line of research, around 1990 
Vilenkin and Shellard proposed that possible effects of cosmic strings in stellar and galactic 
formation and structure should be addressed. 

• Classical/quantum lumps in field theory and elementary particle physics 

The main theme where these highly stable lumps of energy will attract our interest is quantum 
field theory. Many field theoretical models at the heart of our present understanding of ele- 
mentary particles and their interactions have topological defects between the solutions of their 
classical counterparts. Because hadrons, particles interacting via strong subnuclear forces, are of 
two types -heavy (baryons), and light (mesons)- it was tempting to think of them respectively as 
quantum solitons and light quanta. This point of view was pioneered by Skyrme and Finkelstein 
as early as the sixties. The first author even proposed a variation that encompasses solitons 
on the (at that time fashionable) Gell-Mann/Levy sigma model of strong interactions . In the 
Skyrme model, the solitons, usually referred to as Skyrmions, would describe the classical limit 
of baryons whereas mesons were associated with light quanta. 

Needless to say that a puzzling question arose: what is the nature of the quantum field states 
that are the descendants of classical lumps? What do solitons look like in quantum field the- 
ory? The first attempts to explore this territory concentrated on studying the quantum A((/')| 
and sine-Gordon kinks. In 1974 Dashen-Hasslacher-Neveu succeeded in computing the one-loop 
correction to the classical mass of these solitary waves by developing the ?i-expansion of these 
(l-l-l)-dimensional field theories. Moreover, in the second case, where periodic in time soliton- 
antisoliton solutions (breathers) exist, DHN generalized the Bohr-Sommerfeld quantization pro- 
cedure to field theory, obtaining the semi-classical spectrum of these new types of bound states. 
Two years later, Comtet, Cahill, and Glauber provided a closed formula for the expectation value 
of the normal ordered Hamiltonian in quantum soliton states of these one-dimensional systems. 
The CCG formula accounts for the bound states of the second-order fluctuation operator around 
the classical kinks and exactly reproduces the DHN results for static solitons. 

Other techniques for the quantization of non-linear waves were soon developed. To mention but 
a few: 1) Goldstone and Jackiw related the semi-classical expansion to approximations working 
in molecular and many-body physics. 2) Christ and Lee used a collective coordinates method. 
3) Cahill unveiled a variational/coherent state approach. 4) Faddeev and Korepin profited from 
the fact that the sine-Gordon equation is a completely integrable Hamiltonian system with an 
infinite number of degrees of freedom to invent a completely new field: Solving quantum infinite 
systems by means of the quantum spectral transform. 5) Coleman, besides writing a priceless 
review on the subject, showed that the quantum soliton of the sine-Gordon theory was no more 
than the fundamental fermion of the massive Thirring model. Two revolutions were sparked: 
a) Solitons, despite arising in bosonic theories are fermions (like baryons). b) Dualities between 
different models at different regimes of the parameters exist. 6) Mandelstam discovered the 
(non-local) creation operator of the sine-Gordon soliton. 

In the midst of all this excitement, further fuel was added to the fire by three new findings: 
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1) In 1973 Nielsen and Olesen rediscovered Abrikosov magnetic tubes in a different system. The 
Abelian Higgs model supports topological defects that are mathematically identical to Abrikosov 
vortices in a relativistic context. Immediate interest in NO vortices was kindled because they 
were thought of as field theoretical models of dual strings, popular in those days in hadron 
physics. 

2) Looking for a non- Abelian cousin of ANO vortices, also in 1994 't Hooft at CERN and Polyakov 
in Russia independently found extended objects in the Georgi-Glashow model, 't Hooft-Polyakov 
magnetic monopoles are not tubes of magnetic flux but, instead, proper solitons (or point de- 
fects); their energy density is localized mainly in a finite 3D ball with exponentially decaying 
tails, except for an Abelian long-range (-) potential, thus resembling magnetic monopoles from 
afar. 

Abelian ANO vortices have been shown to induce half-integer angular momentum quantum 
numbers on an electrically charged particle and a change in statistics (Wilczek), whereas 'tHP 
magnetic monopoles also carry spin ^ (Jackiw-Rebbi, 't Hooft-Hasenfratz) from a spin/isospin 
mechanism. 

3) In Moscow, in 1975 a Russian quartet -Belavin, Polyakov, Schwartz, and Tyupkin- also 
discovered proper solitons (up to scale invariance) in pure Yang-Mills gauge theory, without 
any interaction with any kind of matter, in (l+4)-dimensions. Because there is no physically 
sensible space-time of 5 dimensions, BPST solitons are considered in 4-dimensional Euclidean 
space. In this context, the fourth coordinate is understood as "imaginary" time, suggesting a 
change of name to BPST instantons and a different physical role: being classical minima of the 
Euclidean Yang-Mills action, instantons dominate the semi-classical expansion of the Euclidean 
YM integral functional. These topological solutions thus provide the leading approximation to 
the tunnel effect amplitude between classical vacua and build the YM vacuum as a Bloch wave. 

Therefore, topological defects dress different physical disguises in different dimensions of the 
space-time in which they live. This, which determines when a given topological defect is a 
domain wall (surface defect), string (line defect), particle (point defect), or texture (instanton), 
lies at the core of the p-brane scan of Townsend. 

• Multicomponent kinks 

Advances in the study of multi-component kinks /solitary waves /domain walls have been achieved 
over the past thirty years. Derrick's theorem forbids the existence of soliton-like solutions in 
scalar field theories with (l+d)-dimensional space-times if d > 1. There are no obstructions, 
however, to the existence of kinks in theories with N interacting scalar fields, provided that the 
space-time is the M^'^ two-dimensional Minkowski space. 

In 1976 Montonen, and independently Sarker,Trullinger, and Bishop proposed a model with 
two real scalar fields and field interactions such that the old kink belongs to the space of 
static solutions of finite energy of this field theoretical model. There was however an important 
novelty: another kink was found, such that the two components of the field profile were not 
zero. To distinguish between the two kinds of solitary waves, the old kinks were denoted as 
TKl - one-component topological - kinks whereas the new kinks were referred to as TK2 - 
two-component topological - kinks. Rapidly, Rajaraman and Weinberg, using the so called 
trial orbit method, identified a special member of a third class of NTK2. The whole manifold 
of non-topological two-component - NTK2 - kinks was identified slightly later by numerical 
integration, but the deep reason for their existence was unveiled by Magyari and Thomas, who 
showed that the system of two ODEs to be solved in the search for kinks is a two-dimensional 
integrable mechanical system: the Garnier system discovered in 1915. The Garnier system is not 
only integrable but Hamilton-Jacobi separable, and Ito took profit from this fact to analytically 
calculate all the kink solutions of the so-called MSTB model. For more than one scalar field, 
simple topological arguments do not ensure lump stability, but Ito and Tasaki classified stable 
and non-stable kinks by using the sophisticated Morse index theorem. Overlooking the difficulty 
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of controlling the spectrum of the second -order fluctuation operator (Hessian), in this case a 
2x2 matrix Schrodinger operator, one of us (JMG) developed the complete Morse theory of 
the configuration space of the MSTB model a la Bott . Models in the same class as the MSTB 
model were addressed by the AAI, MAGL, JMG trio at the turn of the last century. The kink 
varieties were identified and their stability was unveiled in a series of papers. 

(l-l-l)-dimensional models of a complex scalar field with potential energy equal to the square of 
the norm of the gradient of a holomorphic function are very interesting because of the possibility 
of supersymmetric extensions; in fact, these models are obtained by dimensional reduction of 
N = 1 supersymmetric models of one chiral superfield. Vafa et al, in 1989, thoroughly studied 
all the stable kinks arising in these models, whereas Townsend analyzed the balance between 
kink masses. Again, the AAI, MAGL, JMG trio explored the same system by taking a real 
analytic point of view. Another interesting model, in this case coming from the dimensional 
reduction of a AA = 1 supersymmetric theory of two chiral superfields, was addressed by Bazeia, 
Nascimento, Ribeiro, and Toledo -henceforth the BNRT model- in 1995. The kink equations are 
not completely integrable but, Shifman and Voloshin found a complete family of TK2 kink solu- 
tions whereas Bazeia and collaborators studied the stability properties. Although the analogous 
mechanical system is not completely integrable, some of us discovered that for some values of 
the mass of the second boson integrability holds and all the kinks can be found. 

• Recent advances in soliton quantization 

In 1994, the spectacular solution of = 2 supersymmetric Yang-Mills theory by Seiberg-Witten 
provided, as an aside, an exact formula for the quantum mass of BPS monopoles. The same 
result in the low energy domain was derived ten years later, using more down-to-earth meth- 
ods, by the Stony Brook/Wien group formed by Rebhan, van Nieuwenhuizen, and Wimmer. 
This work followed previous investigations by the same team, together with Goldhaber, about 
computations of mass shifts induced by one-loop fluctuations on supersymmetric kinks. The 
extreme elusiveness of this issue did not prevent these authors from identifying the old DHN 
formula as being based on a regularization method that sets a cutoff in the number of fluctu- 
ation modes to be counted, rather than the conventional energy cutoff. Another group from 
Minnesota University addressed the same problem by using high-derivative regularization, with 
SUSY being preserved by boundary conditions to find similar results. Phase-shift analysis by an 
MIT group -Jaffe, Graham, and collaborators- also led to some advances, in this case in a purely 
bosonic setting. Finally, in 2003 van Nieuwenhuizen, Rebhan, and Wimmer, and independently 
Vassilevich, succeeded in computing the one-loop mass shift to the mass of the supersymmetric 
Abelian vortex. 

1.2 A brief history of heat kernel/zeta function regularization methods 

• Zeta function regularization and the heat kernel expansion 

The method of zeta function regularization was invented by Dowker and Critchley and, inde- 
pendently by Hawking, circa 1976. Implementation of standard regularization/renormalization 
procedures in Quantum Field Theory on curved space-time backgrounds led to the introduc- 
tion of this regularization method as the best suited technique to combine second quantization 
phenomena with general relativity. Vacuum expectation values of spatial integrals of the energy- 
momentum tensor are essentially given by the trace of the square root of some differential oper- 
ator of Laplace type. Simili modo, the partition function of Euclidean quantum field theories is 
a functional integral that, up to one- loop order in the ?i-expansion, is the inverse of the square 
root of the determinant of another differential operator of Laplace type times the exponential of 
the Euclidean action over h. 

Traces and determinants of powers of elliptic operators can only be defined by means of a 
process of analytic continuation that mimics the definition of the Riemann zeta function as a 
meromorphic function, giving formal meaning to strictly divergent series in some region (Re s < 
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1) of the complex plane. By replacing natural numbers by eigenvalues (hopefully forming a 
discrete spectrum), generalized zeta functions associated with differential operators are defined. 
There is a general theory of elliptic pseudo-differential operators that characterizes the conditions 
under which the generalized zeta functions are meromorphic functions, and values away from 
poles of the zeta function, and derivatives of zeta, are taken as "regularized" definitions of traces 
and logarithms of determinants of (complex powers of pseudo-)differential operators. 

In interesting physical, cases the pertinent differential operators are those ruling small quantum 
fluctuations in gravitational, Yang-Mills, or solitonic classical backgrounds. Generically, the 
spectral information in these situations is grossly insufficient for identifying the generalized zeta 
function in terms of known spectral functions. Fortunately, B. and C. de Witt had already pro- 
posed, in the mid sixties, use of the high-temperature expansion of the kernel of the generalized 
heat equation provided by the differential operator of Laplace type to unveil the meromorphic 
structure of the generalized zeta function. To achieve this goal, one takes advantage of the link 
between generalized heat and zeta functions via Mellin transforms, such that the residua at the 
poles of the generalized zeta function are proportional to the Seeley coefficients of the heat kernel 
expansion. 

• Heat kernel proof of index theorems 

It is remarkable that almost simultaneously, starting around 1970, parallel ideas were applied by 
Atiyah, Patodi, and Bott to construct the heat kernel proof of the Atiyah-Singer index theorem. 
The context was mathematically much more precise, considering the index of the Dirac operator 
acting on sections of spin bundles tensored with vector bundles on compact spin manifolds 
with or without boundary. The theorem identifies the index of an elliptic operator with some 
characteristic classes of the base manifold: typically the A- genus times the Chern character. 

In contrast to physical situations, generalized zeta functions are well defined because the spec- 
trum of elliptic operators on spaces of sections in bundles with a compact manifold without 
boundary as the base space is discrete. On open spaces, characteristic of physical problems, one 
must impose a rapidly decaying behavior (exponential) at infinity in such a way that the elliptic 
operator will act on spaces of functions. Alternatively, one could consider the same problem 
for manifolds with boundary with spectral boundary conditions a la Atiyah-Patodi-Singer and 
allow the boundary to go to infinity to recover the usual situation in physical problems. 

• Generalized zeta functions and heat equation kernels in physics 

On the physical side, heat kernels and zeta functions have proved to be of great use in the anal- 
ysis of gravitational and gauge anomalies arising in the one-loop approximation to the effective 
action. The computation of quantum effects around gravitational, Yang-Mills, or other classi- 
cal backgrounds using heat kernel/zeta function regularization has been a important theme in 
theoretical physics over the last thirty years. Although many researchers have contributed to 
these developments, we particularly mention the Leipzig/Barcelona group of Bordag, Elizalde, 
Kirsten, Vassilevich and collaborators. 

In particular Bordag and Vassilevich, together with two members of the Stony Brook/Viena 
group, van Niewenhuizen and Goldhaber, used these techniques to compute the one-loop mass 
shift to supersymmetric kinks. Starting almost at the same time, we applied heat kernel/zeta 
function methods to re-work one-loop mass shifts for Xc/)'^ and sine-Gordon kinks in a purely 
bosonic setting. Having established the method, we succeeded in computing mass shifts for 
other kinks in models with a single real scalar kink and non-Posch- Teller Schrodinger operators 
governing quadratic fluctuations. Moreover, the generalization to models with several scalar 
fields having multi-component kinks was reported by us in a series of papers. The generalization 
for dealing with matrix differential operators was the key step that allowed us to compute the 
one-loop mass shift for Abelian self-dual vortices. 
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1.3 Chart of aims 



Research on the quantum descendants of classical topological defects can be classified within two broad 
areas, although with important (l+l)-dimensional exceptions. 

1. In ordinary field theories, the most effective approach is to develop semi-classical analyses or 
^-expansions around the classical soliton solutions, generalizing the old WKB approximation 
method of quantum mechanics to quantum fields. This strategy has so far been fully successful 
in computing one-loop corrections to classical observables only for sine-Gordon and X(p'^ kinks 
and sine-Gordon multi-solitons. 

2. After the seminal paper of Olive and Witten identifying solitons as BPS states in theories 
with extended supersymmetry, taking advantage of this fact much more detailed information on 
quantum corrections to supersymmetric solitons has been acquired. In parallel, the conventional 
semi-classical expansion has been used to estimate the mass shift for SUSY kinks, vortices and 
magnetic monopoles, although great care is needed in combining supersymmetry with suitable 
boundary conditions. 

3. In integrable (l-l-l)-dimensional field theories such as the sine-Gordon system, full information 
on quantum solitons is available due to the existence of an infinite number of conserved charges. 
Also in this case, the identification of solitons as coherent states is enlightening because normal 
ordering is sufficient to achieve full renormalization, and the expectation values of operators in 
coherent states behave as their classical counter-parts. 

Our goal in this set of lectures is to develop semi-classical (weak coupling approximation) analyses 
for multi-component kinks, arising in multi-scalar field theory, and Abrikosov-Nielsen-Olesen vortices, 
arising in the Abelian Higgs model. The one-loop mass shift is essentially the trace of the square root of 
the second-order fiuctuation operator (Hessian), modulo some (infinite) renormalizations. Because the 
spectrum of the Hessian is generally unknown in these cases, we are forced to use asymptotic techniques 
to deal with the generalized zeta function of these second-order matrix differential operators. We shall 
describe our method as applied to multicomponent kinks in Sections §. 5 and 6, whereas a conceptually 
identical but much more technically complex procedure is developed in Sections §. 7 and 8 to compute 
the one loop mass shift of ANO vortices. To explain all the subtleties of our approach in as simple a 
context as possible, in Sections §. 3 and 4 we fully address the problem of computing the (very well 
known) one-loop mass shift of the Xcj)'^ kink. As a bonus, comparison with solidly established results 
obtained by other procedures will provide a precision test for our method. 

In Section §. 2 we offer a summary of heat equation kernels, asymptotic (high-temperature) 
expansions, and generalized zeta functions for a very broad class of differential operators of the type 
that we are going to handle. The connection of these concepts and techniques with the formulas 
arising in our physical calculations is explained in Appendices II, III, and IV. 

1.4 One-loop quantum corrections to soliton masses and the Casimir effect 

The problem of computing quantum corrections to the mass of topological defects is closely related to 
the Casimir effect. The field profile distorts the spectrum of quantum fiuctuations around the ground 
state in a similar manner to the plates of a capacitor in a vacuum. The Casimir effect measures the 
quantum energy of the vacuum when two plates are present with respect to the same quantity without 
plates. The quantum correction to the mass of a topological defect measures the quantum energy of 
the topological defect in its ground state with respect to the quantum vacuum energy. These problems 
lie at the heart of the conceptual foundations of quantum mechanics: there is nothing more quantum 
mechanical than the non-zero energy of nothing!. Throughout these lecture notes we shall refer to 
such things as kink Casimir or vortex Casimir energies by analogy with the quantum energy of the 
Casimir set-up. To justify such an abuse of language, we include Appendix I to describe the Casimir 
effect. 
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1.5 Note on the bibliography 



We shall present the bibliographical References in a global and non-detailed way, except in the cases 
where specific and new results are discussed. It is understood that standard books and monographs 
contain precise bibliographical information. Also, recent References are chosen insofar that they have 
been used in the elaboration of these Lectures. 



• Classical papers, lectures and treatises on solitons 

Important classical papers on the foundations of the matter are: [3], [Ij, and [5]. A complete 
collection of timely mid-seventies works about soliton quantization and semi-classical methods 
can be found in Reference [6] . Seminal lectures on classical lumps and their quantum descendants 
are those of Sidney Coleman in Jaca and Erice, see [7]. In Reference [8] a review is offered 
emphasizing the homotopical nature of topological solitons. The earlier monographic books are 
those of Rajaraman and Drazin [TT]. The Rajaraman treatise aims to address both the 
physically and mathematically relevant aspects of extended states in quantum field theory. The 
Drazin goal is rather mathematical; i.e., the application of techniques of integrable systems as 
the inverse scattering method to find soliton and multi-soliton solutions. An important book 
on vortices and monopoles of the highest mathematical rigor is monograph [TU]. More recent 
treatises such as the books by Vilenkin/Shellard or Manton/Sutcliffe thoroughly address the 
issues, with emphasis on Cosmology in the former, and Mathematical Physics in the latter. 
We also mention the earlier papers on Abrikosov-Nielsen-Olesen vortices |14j and [K] because 
(quantization of) these extended objects is the main concern of these Lectures. 

• Bibliography on generalized zeta functions and heat kernel methods 

The zeta function regularization method started with papers [16] and [T7] in response to the 
need to compute quantum effects on curved backgrounds. Even before it was used as a regulator 
of a physical observable, the generalized zeta function arose as Mellin transforms of heat kernels, 
giving particle propagation in curved spaces. Comparison with particle propagators in Euclidean 
time led B. de Witt, [18], to study the asymptotic high-temperature (short-time) expansion of 
the heat kernel. See also [19] for a modern review on this important mathematical tool. Over 
the last thirty years this broader field, the physical applications of heat kernel expansions and 
generalized zeta functions, has become one of the most important subjects of Mathematical 
Physics. Good References are the monographs [20], [21], and [22]. On the mathematical front, 
we choose |23] and |24j as our favorite References. The link between the coefficients of the high- 
temperature heat kernel expansion and the Korteweg-de Vries conserved charges is explained, 
for example, in |25| . 

• The 1976-1989 period 

This period started with two papers addressing the quantization of supersymmetric kinks |26j . 
[27j whereas the seminal paper of Olive and Witten [28] recognized the link between BPS solitons 
and extended super symmetry. Also, the issue of kink quantization was addressed in the paper 
|3Uj . although in this case it was applied to the exotic Xcj)^ kink discovered in [2H]- Important 
advances in our analytical knowledge of vortex and multi-vortex scalar and vector field profiles 
were achieved in papers [ST], [32], and [33]. In a almost simultaneous development, the MSTB 
model was introduced in papers [M] and |35j . This model is a system of two one-dimensional 
scalar fields having a rich variety of two-component kinks that was first investigated in [36] using 
the trial orbit method. The integrability of the kink equations was unveiled in [37], although 
this fact was not fully exploited until Ito showed that the mechanical system is Hamilton-Jacobi 
separable [38]. The kink stability issue was elucidated in [39] by applying the Morse index 
theorem, and one of us developed the full Morse theory of this problem in [40] and |41j . 

• Recent papers on multi- component kinks 
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Over the last decade many works have been devoted to investigating kink or sohtary wave 
solutions in systems, super symmetric or not, with two or more scalar fields. It is important to 
mention this research because computation of one-loop mass shifts for multi-component kinks was 
the intermediate landmark that allowed us to fulfil the same task for self-dual vortices. Besides 
the MSTB model, which is not discussed in these Lectures, another interesting field theoretical 
model with two real scalar fields was first described in [l2] and [33]. One-component and two- 
component stable kinks with the same energy were discovered very soon after. Shifman and 
Voloshin found in 03] that these kinks belonged to a continuous family, all of them degenerate 
in energy, and hence stable kinks. The AAI, MAGL, JMG trio discovered that, for special values 
of the second boson mass, the static equations are fully integrable and the whole kink variety was 
studied in [50] . The trick is to realize that the mechanical problem is Hamilton- Jacobi separable 
in either Cartesian or parabolic coordinates when the pseudo-Goldstone boson mass is either 
2m or was shown in [35]. Other systems with two scalar fields have been considered, for 

instance in [46] . where kink solutions are discussed in either planar or cylindrical Minkowskian 
space-time. References analyzing kink solutions in models with three scalar fields are [51] and 
|52j . In the case of systems of a complex scalar field, holomorphic superpotentials are naturally 
connected with extended supersymmetry and automatically provide M = 2 BPS kinks, see |47j . 
|48j . and [39], the latter reference offering a thorough analysis of this topic. A recent review 
dealing with these developments and other interesting soliton phenomena is [53]. 

• Recent papers on soliton quantization 

In the second half of the nineties, much attention was drawn to the study of Casimir energies in 
different geometries, see e.g. [53], [55], and [56], a problem close to computing kink ground-state 
energies. The issue of quantum corrections to SUSY kinks was revisited from different viewpoints 
in [H7], [SH], [HS], [SO]) and [HI]. A deeper understanding of the several different regularization 
methods used became available after the work of the Stony Brook/Wien, Minnesota, and M.I.T. 
groups, see also [62] . Another regularization method was applied to the SUSY kink by a Stony 
Brook/Leipzig collaboration based on heat kernel/zeta function methods in [H^]. Almost at 
the same time, several of us applied heat kernel/zeta function technology to calculate one- 
loop mass shifts to the mass of many one-component purely bosonic kinks in the sine-Gordon, 
a variation of the sinh-Gordon, and \4>^ models, see [63]. This work was followed by 
similar calculations applied to one-component and two-component kinks in the MSTB and BNRT 
models in [65], and [66]. After a paper by Bordag on the fermionic vacuum energy in a vortex 
background, [67], the one-loop mass shift to the SUSY vortex was calculated in [63] and [69] . 
In both papers, [70] and [7l], we were able to compute the same quantity for purely bosonic 
self-dual vortices. The one- loop renormalization program in the Abelian Higgs model can be 
found in Reference [72] and is suitable for the goals addressed in this work. In [73] a more or 
less unified formula is offered, giving the one-loop mass shift of kinks and self-dual vortices as 
a truncated series involving the Seeley coefficients starting from the second one. Although the 
Seiberg-Witten solution oiJ\f = 2 SUSY Yang- Mills allows us to know the mass of quantum BPS 
states in any energy regime, the recent interesting papers [75] and [76] provide a more detailed 
knowledge of one-loop mass shifts of M = 2 SUSY monopoles. 

1.6 Note on units and dimensions 

Throughout this work, we shall use a system of units where the speed of light is the unit of velocity: 
c = 1. The Planck constant, however, will be kept explicit because we shall perform semi-classical 
computations. Thus, the dimension of h is [h] = ML, mass x length. These are also the dimensions 
of the Boltzman constant [/c^] = ML, whereas particle masses and temperature have dimensions of 
inverse length: [m] = [T] = L~^. 
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1.7 Brazil lectures 



This work is a written outgrowth of a series of three two hour Lectures given by one of us, J. M. G., 
at the Physics Department of Paraiba University in Joao Pessoa (Brazil) during the third week of 
July 2005. The material presented at each of those Lectures is contained respectively in Sections §. 
3-4, §. 5-6, and §. 7-8 and readers wishing to become acquainted with the physical aspects of semi- 
classical soliton mass shifts can skip reading the rest. We have sketched some brief historical notes 
in the Introduction to place the matter in perspective, at the request of Roberto Menezes, without 
pretensions of completeness or high precision. We also include a Section, §. 2, where heat equations, 
heat kernel expansions, and generalized zeta functions are discussed at a higher level of Mathematical 
rigor. Appendices II, III, and IV are included to establish contact between the spectral functions 
described in Section §. 2 and the physicist's version of the same functions used in the core of the 
Lectures. 



2 Generalized zeta functions and heat equation kernels 

Let us focus on elliptic operators of the general form: 

K = Ko + Q{x) • V + V{x) , Ko = {-A +c^) ■ I , hm V{x) = 

>oo 

acting on the Hilbert space of functions H = ®a=iL\{T'^). Here: (a) T'^ is a toroidal variety, the direct 
product of d circles of radius R= (b) I is the iV x iV unit matrix, (c) V{x) : T'^ Uat^iN) is a 
map from T'^ to the set of iV x iV matrices with real coefficients, (d) Q{x) • V : T*^ ^ T{T'^) (g>Matm{N) 
is a map from T*^ to the tensor product of the tangent space to T*^ times Wlat^{N). (e) V and A = V- V 
are respectively the gradient and Laplacian operators in T"'. (f) is a constant. Assuming that the 
spectrum of K is definite positive, 

Kfnix) = Xn fn{x) , A^, G M > 

the generalized zeta function associated to K is defined as: 

SpccK " 

where s is a complex parameter. Via the Mellin transform 

the generalized zeta function is related to the partition (heat) function h{K) = Tre~^^ of the gener- 
alized heat equation: 



kBT 



The partition function is the integral of the i^-heat equation kernel on the diagonal sub-space of 
T'^ X T'^: 



N N 



Tre"^^ = tr / dvoljd Kk{x,x; P) = [ dvoljd K^^{x,x; p) 

■J A 1 D 1 



A=l B=l 

whereas the kernel itself is the solution of the K-heat equation 

TV 



E (^•<^^^ + i^^^)i^r(x,y;/3)=0^^ , K^^{x,y;0)=5^^.5^'^\y-x) (1) 
with unit source at infinite temperature. 
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2.1 Heat kernel and generalized zeta function for Klein-Gordon operators 

The spectrum of Kq - an N x N diagonal matrix of d-dimensional Klein-Gordon operators- 

B=l 

fc=i 

provides the spectral resolution of the i^o-heat kernel 

K^'^{x,y-P) = 5^^- ^ exp<^-/?( +c2)^-exp^ ^ 



SpecX^ 



AA 



^^-•e-^^^.exp{-^}. exp{-||n(^) + ^|(,^-x)p} 



Spec K^^ 

and the Poisson summation formula 
Y: exp{-t|n(^)+.l2} = (^)i 5; exp|-^^!^^-2vrzK^)..4 ; t = A , v = ^^{y-x) 

leads to the formula: 



/5 



i^^f(x,,1/?)=^^-.e-^^/^.(^)^exp|-^|. exp 
On the other hand, the generalized zeta function is: 

N N d 

yl=l ^(A)g2d ^ R2 ^"^J A=1 fc=l 

Via the Mellin transform, the Epstein zeta function can be written in the form: 

^(^'^'lEl?2e"''^"^) = f7^- E / ^/?/3^-^-exp{-/3( +c^)} 
Again the Poisson summation formula 



E exp{-A.(^(A).^(A))}=^!|!y. J- exp{-^.(/-^^)./^^))} 

allows us to write the Epstein zeta function in terms of the integral representation of Kelvin functions: 

1 /z\-'' 

2 ' \2. 



d ^ 

k=l 

° r('4)eZ''-{o} ° 
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At the infinite volume R ^ oo limit, only the first term survives: 

lim E{s, c'\ f: ^4 ) = vrf iJ'^c^-^^ • ^^^^-1^ . 
k=l ^ ' 

2.2 High-temperature (asymptotic) expansion of the K-heat equation kernel 

In fact, not only when i? ^ oo do Kelvin integrals go to zero but, also, K-y{z) becomes negligible 
when /3 ^ 0, i.e., at the high-temperature limit: 

shows the asymptotic behavior of the "free" heat kernel. 
To find the Ji'-heat kernel, we plug in the ansatz 

TV 

c=i 

in equation ([T|), leading to: 

N 



C=l 

- 5^^-A + g^^(f).9, + y^^(f)}.Cr(x,y;/3) = 0^^ ; (f , y; 0) = (2) 

In the high-temperature limit one can write the heat kernel as the asymptotic series: 

^) •exp{-^-^}.5]5^^5]c^^(f,2/;i^)./3" (3) 

^ ^ C=l n=0 

oo 
n=0 

if C^^{x,y; P) is written as the power expansion above. Plugging ([3]) into ([2]), one obtains the 
recurrence relation between the Seeley densities c^^{x,y; K): 



N r . 

n.6^^ + {xu - yu) ■ - -Qt^m 



E 

C=l 

N 



c^^{x,y;K) 



= i^^'' ■ ^ - Qt''ix)dk - V^^ix)] c^fi(x, y; K) (4) 

starting from: Cq^ {x,y; K) = 6^^. 

Let us introduce the following notation: 

Q01+02H hctd ^^B(;z: -7. 

[X) iim d^<^id^<^2...Q^a, , C„ [X,X,n) 0„ [X) 

d 



< ^ afc < d , ai,a2, • • • ,ad = 0, 1,2, . . . ,d 



k=l 



Thus, in the y ^ x limit the recurrence relations between densities and partial derivatives of densities 
can be written in the compact form: 



15 



k=l 



N ai 02 Qfj / s / 

cti\ / a2 



D=lri=0r2=0 ra=0 



rij \r2 



_ (oi-ri, 02-1-2+1, ■■■,ad-»'d)(^-D-B^^^ + • • • 



+ 



dx^^dx'^ ■■■dx''^ 

N ai — 1 02 



_ (ai-rT_,a2-r2,--,aa-ra+l)(jDB ^^-^ 



+ 



Ql — l\ / 02 



^^2 



ad 
rd 



+5EEE-E''. 

D=l ri=0 r2=0 r£j=0 
gri+r2+ +rd . qAD f^^^ (^ai-l-rua2-r2,-,ad-rd) r^^^ ( 7) -L 
dx\^dxl^---dx'f ^n+lW + 



iv ui U2 — i "d 

4eee-e°= 

D=lri=0r2=0 r£j=0 

9^2^ • • • dx^f 



2 f-—'^ — 1 — 1 — V ^1 / \ ^2 

{ai-ri,a2-l-r2,--,aa-ra)QDB^^^-j _|_ _ 
Af ai a2 ctd-l 

ri,a2-r2,-,ad"l-r-,)(^DB(^) + 



D=lri=0r2=0 r£j=0 

5a;^^ 9x2^ • • • dx^/ 

N ai a2 ctd / 



ad 
rd 



Od-l 

rd 



EEE-E 



D=lri=0r2=0 rtj=0 

to be solved starting from 



riJ \r2 



J dx^^ dx^2 ■ ■ ■ dx^J^ 



_ {ai-ri,a2-T2,---,ad-rd) qDB ^^-^ 



c^^{x,x;K) = 6' 



AB ^ / (-i.«2.-."^)Co^5(x) = 0, if afc / 0,VA; = 1,2, - •• ,d 
(o,o,-,o)(^AA(^) = 1, yl = 1,2, 



2.3 Asymptotics of the partition function and generalized 
zeta function meromorphy 



Defining 



N 
A=l 



the asymptotic expansion of the partition function reads: 



Tre-^^ = e-^'^ 



ttE^ 



n=Q 
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Via the Mellin transform, one writes the generahzed zeta function as the sum of meromorphic and 
entire functions of s: 



Ck{s) = (vri?^). .^trc„(i^) / dp P^^^-'-^ e-'^'^ + —- ■ dpp^-'Tre-^^ 



d \ - 1 



7[s + n-f,c2] 1 



n=0 



r(.) 



+ 



T{s) 



.(5) 



One can show that B{s,K) = ^ ■ dp P'-^ Tre"^^ is a entire function of s (holomorphic in the 
whole complex s-plane C). 



oo , 
d X ^ 1 



n=0 



tr c„(i^) 



7[s + n- ^,c^] 

r(s) 



however, is meromorphic, with poles at the poles of the incomplete Euler Gamma functions: ^[s + n- 



3 The A((/)4)-model on a line 

In the A((/'^)2-model the action 



S 



^29^9^- 4 (2/0,2/)--) 



governs the dynamics of the scalar field ipiuoj v) '■ '^'^ ~^ We choose the metric g^j^y = diag(l, — 1) 
in (l+l)-dimensional M-*^'^ Minkowskian space-time. In our systems of units the dimension of the field 

and the coupling constant are respectively: [V^] = Af 2L2 , [A] = M^^L^^. In terms of non-dimensional 
space-time coordinates and fields 



V2 



m 



m 



the action functional and the field equations of the A((/>)| model read: 



^^2(^0,2;) - ■^{xo,x) = 2(f){xo,x){l - (j)^{xo,x)) 

The shift of the scalar field from the homogeneous stable solution, (j){x^) = 1 -|- H(x^), leads to 
the action 



m 



S = I d X 



2H^{x^') + Ih^x^") 



which shows the spontaneous symmetry breakdown of the internal parity Z2 symmetry. The Feynman 
rules are thus obtained in terms of the Higgs propagator as well as three-valent and four-valent Higgs 
self-coupling vertices: 
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Table 1: Propagator 



Particle 


Field 


Propagator 


Diagram 


Higgs 




iXh 


k 


m?{kQ — hp' — A + ie) 


• • 



Table 2: Third- and fourth-order vertices 



Vertex Weight Vertex Weight 




3.1 Plane waves and vacuum energy 

The general solution of the linearized field equations 



Q^2 (^o,a^) ^^(xo,x) 4(5i/(j;o,a;) = 



governing the small fluctuations of the Higgs field is: 



where fco = = ^/k"^ -|- 4, and the dispersion relation /cq — A;^ — 4 = holds: 

We choose a normalization interval of non-dimensional "length" I = I—tt^, tt^], and we impose 

\/2 2\/2 2v2 

periodic boundary conditions on the plane waves such that = 27rn, n G Z, and the spectral 

density of Kq is: pKQ{k) = ^ = This is tantamount to considering the d = 1, = 1 case of 

Section §.2, although the radius R = ™'h of the spatial circle is slightly modified to fit in with the 

2v 27r 

conventions most frequently used in the literature on kinks. 
From the classical free Hamiltonian 

TTi2\ f r fldSH dSH ld6H dSH , 

k 

one obtains the quantum free Hamiltonian: 
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via canonical quantization: [a{k),a'[q)] = 5kq- The vacuum energy is: 
3.2 The one-loop mass renormalization counter-term 

The Higgs tadpole and the Higgs self-energy are ultraviolet-divergent in the one-loop order of the 
/i-expansion: 

(' (Pk i 

f dk 1 ^2 1 V- 1 /""^^ 

-6i ■ / 4) = -6i • / — • , = ■ -y , = 

y 4vr VFTI mL 2 ^ V . ^ 

A combinatorial factor of ^ has been taken into account in both graphs. The Lagrangian density of 
counter-terms Cc.t. = [<p'^{x) — l) • /(4), giving the vertices in Table 3, must be added to exactly 



Table 3: One-loop counter-terms 


Diagram 


Weight 


X 


6i/(4) 


X — 


6i/(4) 




cancel the divergences above. 
3.3 A(0^)2 kinks 

The configuration space of the classical A((/)^)2-model 

C = {(l){x) E Maps (M,R) < +00} 
is non-connected: C = C+ + U C U C+ _ U C_ + . The energy for time-independent configurations 

^Xj \2dx dx 2^ 

is finite if and only if 



lim = , lim cj}(x) = I , 

rr^itcxD dx x-*±co [ Cp- — — i 

and the four non-connected components of C are classified by the behavior of the scalar field at x = ±00. 
The Bogomolny splitting of the static energy - the energy for time independent field configurations - 



{00) 

i-oo) 
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shows that the absolute minima of E in each sector of C satisfy the first-order equations: 



Besides the homogeneous solutions (j)± = ±1 of energy E{(j)±) = 0, there are kinks or traveling wave 
solutions 



of 




a 



3.4 Kink Casimir energy 



Small kink deformations (t){x) = (pxi-c) +5(j){x) are still solutions of the first-order equations if 5(j){x) G 
KeiD 



Note that 

= D^D 



D6(p{x) 



dx"^ 



d_ 

dx 



+ 2(t)K{x) 54>{x) 



+ 2tanhx ) 6(l)ix) = 
dx I 



+ 4 



cosh X 



K = DD'' 



d^ 



+ 4 



6 



Moreover, the shift of the Higgs field from the stable kink solution, 
to the action: 



(X' 



dx"^ cosh^x 
= (j)K{x) + H{x^'-), leads 



4m^ 



lim 



3\/2A T^ooJ_T 

2 



dxn 



+ 



m 



A 



d^x 



cosh X 



)H\x'') 



2tanhxH^{xf') + -H'^{x''] 



Thus, the Feynman rules in the kink sector must be modified. Both the Higgs propagator and the 
trivalent Higgs self-energy vertex are strongly influenced by the kink. 

The general solution of the linearized field equations 

2 i6H{xo,x)= ^z-^ + K 



— 2-(^o,2;) - -i^{xo,x) + (4 



dx, 







6 



cosh X 



5Xq 



6H{xo,x) = 



for small fluctuations 0(xo,x) = (pxix) + 5H{xq,x) on the kink background 

1 



SH'{xq, x) 



m 



* i\/3: 



v/273 



Ate 



)h{x) + 



m V mL ^ 



1 
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is written in terms of the positive eigenfunctions of the -fT-operator: 



Kfeix) 



dx"^ cosh^x 



k £ 



Eigenvalues 



= 



si = 3 



£2 = A;2 + 4 



Eigenfunctions 



fo{x) 
h{x) 



"T" 



cosh X 
sinhx 



cosh a: 

f^{x) = e^'=^(3tanh2x - 1 - 3iA;tanhx - k'^) 



The choice of periodic boundary conditions in the interval / = [~^^: 2^/2^ 

, mL _ , 
k—= + 5(k) = 2Tm 
V2 



gives the following phase shifts and spectral density: 

, 1 ,mL d6(k) , 



-2arctan 



3k 



2-k^ 



Thus, the sums are over the solutions of the transcendental equations 



71 1 3k 
The classical free Hamiltonian for kink fluctuations becomes 



n e 



2m^ 
hm J 



dx 



d6H d5H 



+ 6H{xo, x)K5H{xo,x) 



dxo dxo 

^ ^V3{AlAs + AsAl) + Y,e{k){A*{k)A{k) + Aik)A*{k)) 

and, after canonical quantization, 

[A3,Al] = l , [A{k),AHq)]=6k, 

one obtains the quantum free Hamiltonian 

Hi2) = (^V3{AIA, + i) + Y,e{k){AHk)A{k) + i)^ 

and the kink Casimir energy 

. ^, , , hm , I- \ ^ hm, ^ ^^1 

AE(0..) = ^(V3 + j:.W) = ^T,./a 



when all the positive modes are non-occupied. 

In sum, the kink semi-classical energy -one-loop order- receives three contributions: 



1. The classical energy, E{(Pk) = 
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2. The kink Casimir energy -zero point energy renormalization- 

AM^ = t.E{4>K) -AEo = ^ f TrK 2 _ TrA'| j 

3. The contribution of Cq.t. to the one-loop kink mass, which is: 

AMI = -3^ • 1(4) . I dx {cPlix) - 4) = 6^ • 1(4) . 

Therefore, the one-loop kink mass shift and the semi-classical kink energy are the divergent quantities: 
AMk = AM^ + AM§ , Es{(I)k) = E{(t)K) + AMk 

4 The kink heat kernel and generaUzed zeta function 
4.1 Zeta function regularization 

We regularize the ultraviolet divergent kink and vacuum energies AE{(j)x), AEq in terms of their 
generalized zeta functions: 

AM^{s) = ^{2^r^,{CUs)-CKo{s)) . 

2 

Here 

and /i is a parameter of dimensions L~^, necessary to keep the dimension of AM^(s) independent 
from the complex variable s. The star means that the zero mode does not enter the kink generalized 
zeta function. Therefore, 



AM^ = lun^ AM^is) = ^ (^(-1) - CKoi-^)) 



and the divergences reappear at s = — ^, which is a pole of AM^(s). AM^(s), however, is a 
meromorphic function of s. 

AM§ can also be regularized in terms of zeta functions. Note that the divergent integral /(4) can 
be expressed as the limit: 



l_ r(s + i) /V^'^^ 

UiL' r(s) '[m^ 



2 

mL mL 



m = - -CKois + i) (7) 



A ..Rr ^ 6/1 /2//2y+2 r(s + i) ^ , 
AM|(s) = -— . . ^ . Cxo(s + 1) 



when the system is defined in the interval [—f^i f^^' "^^^^ 

6h /vy^^ r(s + 
L ' y ■ r(s) 

and 

AM|= lim^AA/|(.s) = ^.Ci.o(^) 
Another, more direct, regularization for 7(4) is possible: 



ri 2//L \ m 



/(4)=lim;^.(-^) .Ci.o(-) ■ (8) 
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The problem is that 



and AM§ = lim^_^i AM§{s) arise at a different point in the complex s-plane from the point where 



AM^ is obtained. 



4.2 The Dashen-Hasslacher-Neveu (DHN) exact formula 

The partition and generalized zeta functions for the vacuum operator Kq are in the R ^ oo limit 0: 



2v/27r 



-4/3 



mL r s--^ ,,_mL 1 r(. - 1) 



The poles of Ci<:o(^) thus the poles of the Euler Gamma function r(s— i): s—^ = 0, —1, —2, • 
The vacuum energy reads: 



AE( 



lim , „ 
2 \ 



22-1 Y{s) 



The partition and generalized zeta functions for the kink operator K can also be given analytically: 

f+oo 1 r+oo 



Tr*e 



-I3K 



-3/3 



mL 

V8^ 



■ e-^^ + e"3^(l - Erfcv^) - Erfc2v^ 



CAs) = CKois) + 



T s) 



+ — / dA: -77-(A;) 



27r 



V^r(s) 



2 ^1 1 3 1, 1 1 



2' 2' 2' 3^ 4«s 



respectively in terms of complementary error functions Erfcx and hyper geometric Gauss functions 
2Fi[a,b,c;d]: 

77 fi I 1 3 1 r(|) ^ i-iy r(/ + i)r(. + ; + ^) 

''^^^2'^^2'2' 3^ r(i)r(s + i) '^^ 3H 



T{1 + 



Thus, besides the poles of Cko{s), Ck{s) has poles at: s + I + ^ = 0, —1, —2, ■ 
The renormalized kink Casimir energy 



-n, • • •. 



hm /2u2V r(e) 
lim 



2^2^ e^o\m'^ J r(-i + e) 



hm 
— ^=— hm 
2V27r '^^o 

hm 



lim 



2V2lT e-*0 



A P fl p 3 _li 3^ 

(-i-|-e)4 2+= 

2;u2 



+ 2 + ln^-31n=^ - 2F[[l,0, f, -|] + o(e) 



2^2 ^ 
- + 3 In — 5" ^ 

e V3 



^In Appendix I it is shown how this limit can be safely taken, leaving no remnants, when PBC are chosen. 



23 



still has a pole; zero-point vacuum energy renormalization is not sufficient. The special values 



1 



vr , 4 
—= — m - 
V3 3 



have been taken into account in the derivation above. 

Does this result agree with the corresponding Dashen-Hasslacher-Neveu (DHN) formula obtained 
via the Stony Brook/Wien mode number regularization method: 



hm 
2V2 



\/3 + — lim 

27r A— »oo 



dk vr 



3hm 

— ^=— hm 
2V2vr A^oo 



A 



dk 



-A VWTi 

The zeta function regularization procedure ([7]) for AMji provides the result: 



6h 



3hm 



lim 



3hm 

— r=— lim 
2^/2^T e^o 

hm 
— ^=— lim 
2^/2vr e^o 



2fi' 



_ /V 

2vr e^'o \ m2 



lim I - 



niL 



1 



8vr 4'^+: 



4-=r(e) 



r(- 



+ e) 



2fi 



- +ln-^ -ln4+ (0(1) 



-+31n^-2(2 + l) 



H-^))+o{£) 



because the difference of the digamma functions is: iJj{1) — 
give the well known answer for the one- loop kink mass shift: 



In 4 — 2. Thus, both methods 



AMk = AMf. + AAf^ 



hm 3hm 



2V6 ttV2 



In the DHN derivation, however, when the mode number regularization cutoff is used, the second 
(negative) summand in the formula comes from the kink Casimir energy AM^. The first summand is 
due to the non-exact cancelation between the ultraviolet divergences arising in AM^ and the induced 
energy AM§ by mass renormalization added to the contribution of the bound state. In our zeta 
function regularization procedure, the origin of the two terms is more clear: the first summand comes 
from the finite piece in AM^ at the physical point s = —■^ but the other piece is found in the regular- 
ization of AM^; i.e., AM^ does not exactly cancel the divergence of AM^ but the renormalization 
process leaves finite reminders in the kink Casimir energy and the mass renormalization counter-terms 
induced, providing the correct answer. 

The alternative regularization of 7(4) ([8]) applied to AM^ leads to the result: 



am| 



^ lim ( 

3hm 

, lim 

ySvf £^0 

3hm 

hm 



2^/2vr e^o 

hm 
— ^— lim 

2V2vr e^o 



mL 



1 



8vr 4" 



2/i^V 4-^r(g) 

r(i + e^ 
2/x2 



Tis 



Tis) 



+ In 



m^ 



ln4+ (-0(1) 



+ 31n 



2n^ 



m^ 
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because the difference of the digamma functions is: V'(l) ~ ^(i) — l'^'^- Thus, using the latter 
regularization for 7(4) we would obtain 

\Mk = AM^ + AM| = 

This (bad) result was achieved in the literature on the matter by regularizing the ultraviolet divergences 
by means of a cutoff in the energy, rather than in the number of modes. Also, one could give this 
answer without taking into account the zero mode in the CCG formula, see section §. 6.2 . 

4.3 The high-temperature expansion of the partition function 

Even without complete knowledge of the spectral data of the kink fluctuation operator it would 
be possible to obtain an approximate formula for the one-loop mass shift from the high-temperature 
expansion of the partition function. 

The heat equation kernel for the KQ-heaX, equation 
d d'^ \ 

^ - ^ + 4 1 K/ro(^>y;/5) = , KKa{x,y;G) = 5{x-y) 
of the vacuum fluctuation operator Kq = — ^ + 4, for (3 small, is: 

The kink fluctuation operator K is the Schrodinger operator 

cP 6 
K = —-. + A + V{x) , y(x) = 60|(x)-6 



whereas the corresponding K-heaX. equation kernel 

(^-^+4 + y(x))/^i^(x,y;/?)=0 , KK{x,y;Q) = 5{x - y) 
can be written in the form 

KKix,y;(3) = KKo{x,y;f3) • CK{x,y;l3) 
if CK{x,y; P) satisfies the transfer equation 
d X — y d d'^ 



d(3 f3 dx dx^ 



+ V{x)]CK{x,y;P)=0 , y; 0) = 1 



and it is set to be unity at infinite temperature. 

Solving the transfer equation as a power series in /? 

oo 

CK{x,y;P) =^ Cn{x,y;K) , co{x,y;K) = l 

n=0 

the PDE equation becomes tantamount to the recurrence relations: 

(n + l)cn+i{x,y;K) + {x - y)^^^^{x,y; K) + V{x)cnix,y; K) = ^-^{x,y;K) 

ox ax^ 
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The high-temperature heat equation kernel for K is thus given as: 



KK{x,y;P) 



-4/3 



n=0 



We are actuaUy interested in the trace of the heat kernel to find the partition function for small /?. 
The recurrence relations become 



Cn+i{x,x;K) 



1 



n + 1 



^^^Cn{x)-V{x)cn{x,x;K) 



when limj^^^:- To deal with this delicate limit, we have introduced the following notation: ^^^Cn{x) = 
. RecaU that ^'^)Co(x) = lim^^a: = 5''^. We also need (obtained after differenti- 



lim. 



ating the first recurrence formula /c-times) recurrence relations among derivatives: 

k 



1 



n + k 



The high-temperature asymptotic expansion of the partition function reads: 

^-4/3 oo 



n=0 



mL 
2\/2 



Cn{K) = lim / dxCn{x,x;K) 

L— >00 J mL 
2V2 



Using the recurrence relations the c„(x, x; K) densities can be found -they are the conserved charges of 
the KdV equation, see Appendix II- and, via integration over the whole line, the kink Seeley coefficients 
are obtained: 

'°<^'> = £"„7f • -"(^'= (2n-l)!! ■ ■ 
4.4 The Mellin transform of the asymptotic expansion 

To obtain the generalized zeta function from the asymptotic expansion of the partition function, the 
Mellin transform is split into two integrals, inside and outside the convergence radius: 



Ck{s) 



1 



r(.) 



1 00 i-l poo pi 

^= ■ y Cn{K) ■ / d(3 2 . e-^^ + / dl3 13'^^l\*e-^^ - / 
v4vr ^ Jo Ji Jo 



s-l 



On general grounds, it is possible to show that 

1 

Bk{s) = -— ci/?/?-iTr*e-^^ , 

r(s) Ji 

where the star means that the zero eigenvalue is not accounted for, is an entire function of s. 

Note, however, that the zero mode is included in the heat kernel high-temperature expansion. 
Subtraction of the contribution of the zero mode to C^(s) coming from the high-temperature range of 
the Mellin transform is a tricky affair. In fact, 

^1 

d(3(3'-^ 



is an improper integral if Res < 0. We define this integral in the spirit of zeta function regularization 
as: 



lim -1- . r dl3l3'-^e-'f' = lim 1 • #4 
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Near the physical value e = 0, the behavior of the Euler incomplete Gamma function is such that the 
regularized integral is: 



1 



1 



s+l 



1 



1 



sT{s) T{s + 1) 



The zeta function regularization procedure directly provides the value for the improper integral that 
would be obtained if the divergent integral I had been renormalized by adding another (related) 
divergent integral: 

1 



tR 



, X lini 
Vis) c-^O 



A similar strategy must be adopted for computing the zeta function of the vacuum (Klein-Gordon) 
operator to find: 



Ci^o('S) 



rah 1 



r(s) 

mL 1 



Svr r(s) 4^-^ 



7[s-i 4]+r[.-i 4] 



The incomplete Euler Gamma function j[s — ^,4] has poles at s — ^ = 0, — 1, — 2, — 3, • • • but its 
complementary r[s — ^,4] is an entire function of s. 

By the same token the generalized zeta function of the kink fiuctuation operator reads: 



Ck{s) 



T{s) 



No 



7[s + n-i4] , 1 



n=0 



+ 77^- E W 



7[s + n-i4] 1 



n.=Afo+l 



+Bk{s) 



7[s + n — ^,4] has poles at s + n — ^ = 0, —1, —2, —3, • • • and a large but finite number A'^o is chosen 
to separate the contribution of the high-order coefficients. 



n=No+l 



is holomorphic, however, for Res > —Nq — 1. 



4.5 The high-temperature one-loop kink mass shift formula 

Neglecting the (very small) contribution of the entire functions, the kink Casimir energy becomes 

AMS^ ~ - • lim f • ^ • T^T- 

^ 2 s-.~\\m^ ) ^ T{s 

i.e. the zero-point vacuum energy renormalization takes care of the term coming from cq{K). 

The other correction due to the mass renormalization counter-terms can also be arranged into 
meromorphic and entire parts: 

7[s + ^,4]+r[s + ^,4] 

The mass renormalization term exactly cancels the ci{K) contribution. Our minimal subtraction 
scheme fits the following renormalization prescription: for theories with only massive fiuctuations, the 
quantum corrections should vanish at the limit where all the masses go to infinity. 



No 



7[s + n- i,4] 



n=l 



\s+n- 
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We end with the high-temperature one-loop kink mass shift formula: 



hm 



^ 2^ ^"^^^ 4n-l +^ 



n=2 



A good test for the appropriate regularization prescription chosen for the zero mode subtraction is 
the value 

(0) _ hm 



AM 



K 



2v/27r 



-0.199471^m 



in perfect agreement (up to a factor due to different conventions) with the result obtained by 
Glauber, Comtet, and Cahill, in 1976. 



4.6 Mathematica calculations 

Computational limitations impose a practical bound on the choice of Nq. Knowledge of, say, ^^^C2 
requires computation of 9 densities: 

(4) Co Co (2) Co (i)Co (o)Co 
(2)Ci (i)Ci WCi 

(0)C2 

In general, the evaluation of ^'^^Cn{x) requires previous calculation of 



1 + 3 + 5 + 7 + • • • + 2?i - 1 + 2?i + 1 



in + 1) 



('=)C7o(x) densities. In fact, because (o)Co = 1 and (^^Co = (2)^0 = • • • = (2«)Co = there would be a 
need to compute only coefficients, but the computer ignores this circumstance. 
Observe that 



AMk = -0.199471/im + DN^hm 



D 



No 



n=2 



8v/27r4 



n-l 



is far from the exact result without adding DN^hm: AMk = — 0.471113?im. 

In the following Table Mathematica, calculations of the Seeley coefficients and partial sums D^^ 
are shown up to A^'o = 10 



n 


Cn{K) 


No 




2 


24.0000 


2 


-0.165717 


3 


35.2000 


3 


-0.221946 


4 


39.3143 


4 


-0.248281 


5 


34.7429 


5 


-0.261260 


6 


25.2306 


6 


-0.267436 


7 


15.5208 


7 


-0.270186 


8 


8.27702 


8 


-0.271317 


9 


3.89498 


9 


-0.271748 


10 


1.63998 


10 


-0.271900 



giving the very good result: Diq = — 0.271 900^m. Finally, we find 

AMk = -0A7137lhm 
with an error with respect to the DHN result of 0.0002580/im. 
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dp[ + 3 + 3 — + 



2^2^ Ji \ 2/3i 2/?t 2/3t V^^P J 

« 0.00032792nm 

is almost the total error. The deviation from the total error is — 

^6^0 (_i) ^ lo^^hm. 



5 The BNRT-model on a line 



In this Section we analyze a model in (H-l)-dimensional M^'^ space-time of two scalar fields with 
dynamics governed by the action functional: 



S 



This model was originally discussed by Bazeia, Nascimento, Ribeiro, and Toledo (BNRT). The main 
novelty with respect to the model studied in the previous Section is that there are two real scalar 
fields in this system that can be assembled into a "isospin" vector field: 

2 

^(yo, y) = X] ^'^(2/0, y)ea ■ K^'^ ^ , e; • 4 = (^at , a, 6 = i, 2 , 

a=l 

where ei , 6*2 form a orthonormal basis in the target space M? . The dimensions of the fields and the 
coupling constants are respectively: [t/ja] = M^L^ , [A] = [u'^] = M^^L^"^, and [m] = L^^. In terms of 
non-dimensional space-time coordinates, fields, and parameters 

m VA A 

the action functional and the field equations of the BNRT model read: 



4m? 



j 1 ^ dcl)a dcl)a 1/. ,2 , o ,2 
V. a=l ^ 



-(xo,x) - -^-^{xq,x) = 2(/)i(xo,x)(l - 2a{a + 1)4)1{xq,x) - A(t){{xQ,x)) 



dx"^ ' dx"^ 



-^^(xo,x) - -^^(xo,x) = (T(/>2(xo,x)(l - 2(T(/>2(xo,x) - 4(cr + (xq, x) ) 

There are four homogeneous stable solutions: 

1. (^±«(xo,x) = ±^-ei 2. 1 .g-.^ 

^ V 2(7 



To quantize, we choose the ^^^"^^ vacuum (because it is the asymptotic value of the generic kinks of 
the model) and shift the fields from it 

</)i(x'^) = ^ + i7(x'^) , i>-,{x^) = G{x^) 
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to write the action in the form: 
.2 



Atd ^ C 

- (fx {m^ixf") + 2a{a + 1)H{x'')G\x'')} - 

_ ^ J j^2H\x^') + 2(7{a + l)H\xf')G\x^') + ^G\x^')^ 

showing the spontaneous symmetry breaking of the Z2 x Z2 vierergroup generated by the internal 
reflections — > —(pi and (p2 — > —4'2- The Feynman rules are thus obtained in terms of Higgs and 
(pseudo) Goldstone propagators and three- valent and four-valent vertices: 



Table 4: Propagators 



Particle Field Propagator Diagram 



Goldstone G{x) 



Am?{kl - /c2 - 4 + ie) 
i\h 

Am?'{kQ — A;2 — + ie) 



Table 5: Third- and fourth-order vertices 



Vertex 



Weight 



Vertex 



Weight 




-96i- 




-192i- 



m 
liX 



Tfl 
An 



-32a{a + l)i 



m 



Xh 



5.1 Plane waves and vacuum energy 

The general solution of the linearized field equations 



— -^(xo,a;) —y-[xo,x) + 4dH[xo,x) = 

OXq ox'^ 

d^6G. , d^5G. , 2.^, N n 
-^^{xo,x) - -Q^i^o, X) +a dG{xo, x) = 
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governing the small fluctuations of the Higgs and (pseudo)Goldstone fields is: 



onixQ,x) = -A/ — ^ , 



5G{xo,x) 



a/a / h \ - 1 r / N 



where ko = uj{k) = Vfc^ + 4, go = 7(5) = y^g^ + cr^, and the dispersion relations /cg ~ ~ 4 = 0, 
- g2 _ = hold: 



„ikx 



1 ^0 \ giqx 



r,igx 



dx'^ 



+ 4 







We choose a normalization interval of non-dimensional "length" mL, I = [—^^, ^^], and impose PBC 
on the plane waves so that: k ■ mL = 27m , q ■ mL = Inr with n,r G Z. Thus, the spectral density of 
Kq is: 



dn r\ 

" dk 



1 



-mL 



vr 



This is tantamount to considering the d = \, N = 2 case of Section §.2. 
Prom the classical free Hamiltonian 



\ J 1 2 V dxf) dxQ 2 dx dx J 2 \ dxo dxo 2 dx dx J 
2 

+ 26H ■ 6H + y • 6G 

= Y.^"^ [^{k){a\{k)ai{k) + ai{k)a\{k)) + ^{k){al{k)a2{k) + a2(A:)a^(A:))] 
One goes to the quantum free Hamiltonian: 

via canonical quantization: [a;,(/c), al:(g)] = 5bcSkq- The vacuum energy is: 
5.2 One-loop mass renormalization counter-terms 

There are three ultraviolet divergent graphs in one-loop order of the /i-expansion contributing to: 
• The Higgs boson tadpole: 



-12i ■ 1(4) - 2cr{a + l)i ■ I{a'^) 



-12i • 



£k 



(27r)2 (A;2 _ fc2 - 4 + ie) 



2a{a + l)i 




+ 



d'k 



(27r)2 (fc2 - A:2 - ^2 + ie) 
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• The Higgs boson self-energy 



-24i ■ 7(4) - 4o-(a + 1) • /(cj^ 




. f dk 1 ^ , , f dk 1 

The Goldstone boson self-energy: 



47r x/FTa2 




-Aa{(j + l)i • 1(4) - 6a^i • /(a^) = 

M/" dA: 1 n f dk 1 

' J 47r ^/WTl J 47r ^/WTo^ 

Due care is needed to take into account a combinatorial factor of ^ in all these graphs. The Lagrangian 
density of counter-terms: 

Cc.T. = ^[Q-m + a{cj + l).I{a^)\[Acfi{x^)-l) + 

+ ^[2(a + l)./(4) + 3a./(a2)]2a<^i(x'^) , 

giving the vertices in the next Table, must be added to exactly cancel the divergences above. 

Table 6: One-loop counter-terms 



Diagram Weight 

X 2i(6/(4) + a{a + l)I{<y')) 

)^ 4i(6/(4) +a(<T + 1)7(^2)) 

-)<J 2ia{2{cj + 1)1(4) + 3a/((T2)) 



5.3 Moduli space of BNRT kinks 

The configuration space of the classical BNRT model 

C = {(t)a{x) G Maps(R,M2)/^[(/,^] < +00} 

is non-connected: 

c = c\:^^vAC[_^_uc\.'^_vAC[_-^^ , I, J = 1,2 

The energy for time-independent configurations 
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is finite if and only if 



lim ^ = 



lim <p{x) 

a;— >itoo 



and the sixteen non-connected components of C are classified by the behavior of the scalar field at 
X — > ±00. 

The existence of the superpotential 




W((^) 




Figure 1: The —U{(j)i, 4>2) potential (left) and the superpotential W{(j)i,4>2) (right) 



allows for the Bogomolny splitting of the static energy: 
2 



a=l 



^ 4m3 /■ , 1 



dW 



\ dx d(j)a J V dx d(f), 



dW 



{W{ct>a{+^))-W{(t>a{-^))} 



showing that the absolute minima of E in each sector of C satisfy the first-order equations: 



dx 



•a ^dW_ 



, a, /3 = 0, 1 



(9) 



(10) 



5.3.1 Kink flow lines 



The solutions of (llOp are the flow lines of the gradient of W ^ and those starting and ending at critical 
points of W are the kink orbits. 










*' 1 












Figure 2: The -U{(f>) 



'^^^^^ ~ ^f^f^ potential (left). Flow-lines: in the ranges c G {—oo,c^) (middle 



left), c = c (middle right), and c € (c , 00) (right) 
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The kink flow lines can be obtained analytically by integrating 

2 

to find, after use of the integrating factor |(/'2|~°^, 

,9 cr ,9! c,,,2 1+'^ 



^+ln|02|)=^ ,1=1, cG (-00,-1 + In 2) 



Thus, there exits a family of two-component topological kinks -TK2 kinks- parametrized by the inte- 
gration constant c, all of them with the same energy: 

Note that there is a critical value = | j-^(2o")~^ if cj 7^ 1, or, c*^ = — 1 -|- ln2, if o" = 1 beyond 

which the flow lines go to infinity and do not correspond to kink orbits because Vl^(|(/)|(oo)) = 00 and 
the energy becomes infinite. 

In general, the kink profile, the dependence on x of the fields for a given kink orbit, cannot be 
expressed analytically. There are, however, two special kink orbits for which this is possible: 



b™{x) = (-l)"itanh(x - a) 
(g^^{x) = 



00 , S i,TKwiN n ; a = 0, l,aG 



J 0TK2(^) ^ itanh[2(l -a){x- a)] 

' \ <^r^(:.) = (-ir/i^sech[2(l-a)(.-a)] ' "i'"2 = 0,l • 

If c = —00, the orbit is a segment on the abscissa axis and one finds the very well known A((/))2 kink 
buried in the kink variety of the BNRT model. Since only one of the two components of the scalar 
field is different from zero, this kind of kink is termed TKl kinks -one-component topological kinks- 
in the context of the BNRT model. 

When c = things become more interesting. The kink orbit is a half-ellipse (if a 7^ 1) and the 
two components of the scalar field are not zero for this kind of topological -TK2- kink. Note that in 
both cases the kink profile is a function of a real parameter a G M, obeying the freedom of setting the 
kink center. This is a general feature shared with other kinks for which analytical expressions giving 
the kink profiles are not available. Therefore, the moduli space of kinks in the BNRT model is the 
two-dimensional manifold M.\ = (— oo,c'^] x (—00,00) with coordinates (c, a). That is, each kink in 
the BNRT model is determined by its orbit (c G (— oo,c'^]) and its center (a G (—00, 00)). 

We close this subsection by enumerating the main features of non-generic kink orbits. 

1. c = — 00: The second component of the scalar field is zero and the kink orbit is invariant under 
the 02 —(t>2 transformation. These orbits belongs to the C}^ or sectors. 

2. c = 0: The kink orbits are half-ellipses also living in the C}^ or sectors. 

3. c = C': These orbits are the boundary between bounded and unbounded fiow lines. They belong 
to or C|i sectors. 
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5.3.2 Kink profiles from integrable systems 



The search for finite energy static solutions in (l+l)-dimensional theories of two scalar fields is tanta- 
mount to investigating finite action trajectories in a mechanical system of two degrees of freedom. The 
potential energy of the mechanical system is equal to minus the potential energy density of the field 
theory. The mechanical analogue in this sense to the BNRT model is a Liouville integrable system 
if fj = ^ and a = ^. Liouville systems are Hamilton-Jacobi separable and all the trajectories can be 
found. Here we shall describe only the cr = ^ case, where the HJ equation is separable using elliptic 
coordinates. Standard application of the HJ procedure provides the following formulas for all the kink 
profile^: 



-1) 



1 sinh ((x — a)) 

2 cosh {{x — a)) + 6^ 



0™[x;a,6] = (-1) 



Ct2 



^JlP' + cosh ((x — a)) / 

The variety of kinks depends on two integration constants: a,b £ (—00,00). The second one, 5 = 
±yj is determined by the constant c G {—00,^), giving the kink orbit. In Figure 3 it is observed 

that upon increasing b a splitting into two kinks arises. 




Figure 3: Solitary waves corresponding to: (a) 6 = 0, (h) b ~ VO.5, (c) b = 1 and (d) b = v30. 

This phenomenon is better understood by studying how the energy density varies with b. The 
energy density 



c-TKr c^u-S" _ ^ + 76^cosh[x] + 46^cosh[23;] + 62cosh[3a;] 

[x;0,6J -2_^^^.^^ - 2(62 + cosh[x])4 

is shown in Figure 4. 



a=l 



Figure 4: Energy density 0, b] for (a) b^O, (b)b^ V0.5, (c) b ^ 1 and (d) b ■ 



30. 



The critical points of S'^^lx; 0, b] 



ds 



TK 



[x;0,6] 



2 sinh X 



P3(coshx) 



d8- 



TK 



■[x;0,6] = 



dx ^ ' ' -)- coshx)^ ^ ' dx 

are the origin x = ~ sinhx = 0, Vb^ and the roots of the third-order polynomial P3(coshx): 



P3(coshx) 



-6^ cosh^ x-b" cosh^ x + (-36^ + 46^) cosh x + hb" -4. 



^There is an open problem related with the search of kink profiles in this model: for (t = 3, (7=i,a = 4, and cr = \ 
the kink profiles can be analytically expressed in terms of elliptic functions, although these solutions are not given in the 
literature. 
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which under the change of variables P3(coshx) = —b'^P{u), v?{x) = — 1 + coshx becomes: 
P{u) = {u^ + {h" + 3)(n2)2 - (46^ - _ Q){u^) _ 4 (^6^ ^ 6^ _ ^ _ JL 

In sum, use of the Cardano-Vieta formulas to solve cubic algebraic equations prompts us to the 
following conclusions: 

1. If 6^ € [0, 1] there are no real roots P3. The only critical point of the energy density , a maximum, 
is the origin and only one lump of energy is carried by these kinks: 

2 



(62 + 1) 



2. If 6^ > 1^ p{u) = has two real solutions: u± = ±^Jr{b'^) £ M. The points 

X = ±771(52) ^ m{b^) = ^— ^arccosh(l + r{b'^)) 

are the maxima of S"^^ on the real line (the origin is now a minimum) and these kink profiles 
carry two lumps of energy . Understanding these two lumps as fundamental particles, a is the 
center of mass whereas b is the relative coordinate of this system of two particles that become 
glued when b^ < 1. 

5.3.3 Kink profiles: numerical methods 

For generic values of a one must rely on numerical integration methods to find the profiles of the 
kink solutions. We do this by solving the first-order equations by standard numerical methods with 
"initial" conditions: 

^i(o) = o , ^(;^/2i0)-^\M0)\i = \ . 

The reasons for this choice are twofold: (1) For any kink solution, 4>i{x) always has a zero. Transla- 
tional invariance allows us to set the zero at x = 0; (2) To ensure that we will find a numerical kink 
solution, we fix 02(0) on a kink orbit for a given value of a and arbitrary choices of c. 

The numerical method provides us with a succession of points of the kink solution generated 
by an interpolation polynomial. The plots of the numerical results show that the behavior derived 
analytically for the kink profiles when cr = ^ is generic. For any value of a the kink profiles are 
composed of two kinks. The parameter c giving the orbit is related to the kink separation. In some 
range of c the two kinks melt into a single kink. The precise value of c at which this happens depends 
on the value of a. a = ^ is singled out, because in this case c = is the value where two kinks fuse 
into a single kink, or, viceversa, a single kink splits into two kinks. 



5.4 TK2 kink Casimir energy 

Small kink deformations (pai^) = (p^^[x;c) + d(j)a{x) are still solutions of the first-order equations if 
6(j)a{x) e Ker Z)(c): 

-£-AM^;c) -2aM^;c) \ ( SM^) \ ^ ( \ 
-2aM^-c) -j^-2a4>i{x;c) ) \H2{x) ) \Q) 

We shall use the notation (p^^{x;c) = (j)J^{x;c)ei + (f>2^{x;c)e2 because no analytic expressions for 
the kink profiles are known (except for some special values of a). The c parameter tells us what kink 
orbit is chosen. Note that: 

+ 84>l{x;c) +Aa{a -l)4>l{x;c) + 2 8a4>i{x; c)Mx] c) \ 

8aMx;c)(f2ix;c) + 4a{a - l)^l{x;c) + 2a^^l{x;c) + a j 



D{c)5(t){x) 
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K{c) 



+ 2#^(x; c) + iaia + l)</>^(x; c) - 2 



8cr(fj + l)0i(x;c)(/'2(a;;c) + 4o-(o- + (j;; c) + 6o-^(/)^(x; c) - cj 

if a = /3. For a / /?, /^^(c) = L»t(c)D(c) and K{c) = D{c)D\c) are exchanged. 



8(t((T + l)</'i(2;; c)(/)2(x; c) 



Moreover, the shift of the Higgs and Goldstone fields from the stable kink solution, <j){x^) = 
{x; c) + H[x^)ei + G{x^)e2, causes the action in the kink sector to be the complicated expression: 



4m^ 
= — lim 

3A T^oo 

+ — I d^x 



^ 2 

axo + — 

T A 

2 



cPx 



-df^Hd^H - (l2(/)f (x; c) + 2(j(a + 1)02(2;; c) - l) H^ix^") 



1 



^d^Gd^'G - ( 2a{a + l)(/>^(x; c) + 3a^^l{x; c) - ^ ) G^(x' 



d^x [8a(CT + l)(/)i(x; c)(l)2{x; c)F(x^)G(x^) + 



m 

T 



(f^x [8(^1 (x; c)H^{x^') + 4a((j + l)0i(x; c)/f(x'')G2(x^)+ 

+4a{a + 1)02(2;; c)f2(x'^)G(x^) + 2a^Mx; c)G''^{x'') 



Both the Higgs and Goldstone propagators, as well as all the trivalent vertices, are distorted by the 
kink. In the background of a TK2 kink, Higgs and Goldstone particles can transform into each other 
by means of the bivalent vertex induced. 

Thus, the classical energy for small fluctuations 0(xo,x) = i^^(x;c) + 6H{xo,x)ei + 6G{xo,x)e2 
reads: 



^(2) 



dx 



d6H dSH 



dxo dx 







+ 6H{xo,x)Kii{c)6H{xo,x) + 6H{xo,x)Ki2{c)5G{xo,x) + 



SG{xo,x)K2i{c)6H{xo, x) + 6G{xo, x)K22{c)6G{xo,x) 



where Kab{c) are the matrix elements of the second order fluctuation operator, which we rewrite, 
together with his supersymmetric partner, in the form: 



K{c) = D{c)D\c) = 
K~ = D\c)D{c) 



-^ + 4 + Vii(x;c) 



^12(2;; c) 



^21(2;; c) 



+ + 1/22(2;; C) 



dx 



+ 4 + T/„(x;c) 



dx2 ' ^ ' •^llV'^l' 

V^^{x;c) 



1/12 (2;; c) 
-£i + cy^ + V^^{x;c) 



240?(x; c) + Aa{(T + l)0|(x; c) - 6 
8(T(cr + 1)01 (x; c)02 (2;; c) = V21 (2;; c) 
4cr(cj + l)0?(x; c) + 60-^02(2;; c) - a{a + 1) 
80? (x; c) + Aa{a - l)$i{x\ c) - 2 
8(T0(x; c)02(x; c) =y27(2;;c) 
4cr(cr - l)0?(x; c) + 2cr202(x; c) - a{a + 1) 

For instance, the family of Hessian operators for the topological kinks found in the cf = \ case can be 
written explicitly: 



Vii(x;c) 
1^2(2;; c) 
V'22(2;;c) 
V\7(x;c) 
V\2(2;;c) 
1/22 (2;; c) 



K{b) 



+ 



3b 



+ 



6 sinh X 



cosh x+b (cosh x+b)'^ 
G-^/fesinh X 
(cosh x+b) 5 



d2 



(coshx+fe)^ 



dx'' ~^ 2 cosh x+b ~^ 4 (coshx+6)-^ 



+ 



3 sinh X 
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Figure 5: Diagonal components ^11(2;) (red) and ^22(2;) (green) of the potential for c=-30, c~-l, c~0.1, 
0=0.245 and 0=0.249. 



In any case, the spectral resolution of 

K{c)fe{x)=e^Mx) , lim K{c) = Ko 

x—>±oo 

for any value of a has the following features: 

1. The kernel of K{c) is of dimension two and the orthogonal basis is provided by the two eigen- 
functions (zero modes): 

/ 90^/ ^ \ / 5<^™^ ^ 

obeying the fact that motion in the kink moduli space costs no energy: all the kink solutions 
are in neutral equilibrium. 

2. There are bound states in a number, N{c), that depends on c. 

/r(^) 

3. There are two branches of scattering states: 
(1) e^ = fe^ + 4 ; #(x;^))=e--f , lim f ^^^^'^^^'f ^ 

.±1 



ez>0, / = l,2,...,iV(c) ; f,,{x) = { ] , I dx f,\{x)fs^{x) = 6'' 



In terms of non-explicitly known phase shifts -determined by (k)-, Periodic Boundary Condi- 

rnL mL 
2 ' 2 



tions in the interval / = [— ^] give the spectral densities in Spec/C(c): 



k ■ ml + 6i(k) = 27rn , pl^fA:) = — (mL + ^f^) 

ZTT dk 



The phase shifts are now read from u^{q) and the PBC give the spectral densities in the other 
branch of Specif (c): 

q-mL + 62{q)=27rn , p'^^(q) = l.{mL + . 
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The eigenfunctions in the continuous spectrum of K{c) also satisfy orthogonahty conditions: H 

dx/tW(x;A:)/y)(x;g) = 5"5fc, , I, J =1,2 . 
The general solution to the linearized field equations 

{2A^i{x; c) + 4a{a + l)^l{x; c) - 2)6H{x^') + Sa{a + l)4>i{x] c)^2{x; c)SG{x^') = 
+ 8a{a + l)Mx; c)Mx; c)6H{x'') + (4(o-(o- + l)4>lix; c) + Qa'^4>l{x; c) - a)6G{x^') = 
is a linear combination of the eigen-functions of K{c) 0: 



d'^6H 


d'^dH 


dxl 


dx"^ 




d^5G 


dxl 


dx"^ 



5H'{xq, x) 
6G'{xo, x) 



+ "^-jl=^A^^\k)e-'^''''^''''fP{x)+A<^\k)e''^^^^^ + 



9 

The classical free Hamiltonian for kink fluctuations becomes: 



^ I ^^"^ 



and, after canonical quantization, 

one obtains the quantum free Hamiltonian 



and the kink Casimir energy 



AEi^^ic)) = ^ ( E + E ^(^) + E ^(«) ) = ^Tri^(-)^ 
\i=i k 1 / 

when all the positive modes are non-occupied. 

In sum, the TK2 kink semi-classical energy -one-loop order- receives three contributions: 



^There is a very subtle possibility. If £jv(c) = sik = 0) or \(q — 0), i.e. if the last eigenvalue in the discrete spectrum 
coincides with the threshold of any branch of the continuous spectrum, half-zero modes enter the game. The Levinson 
theorem in one dimension forces us to include a weight of ^ in the contribution to the energy of those states. 

*The eigenfunctions belonging to the kernel of K{c) are excluded because they do no contribute to the energy. The 
prime in the fields refers to this exclusion. 
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1. The classical energy E{$^^{c)) = 

2. The TK2 kink Casimir energy -zero point energy renormalization- 

AM^(c) = A^((^K(c)) - A^(0±(i)) = ^ (^K{c)^2 - TVk| 

3. The contribution of Cct. to one-loop TK2 kink masses is: 

AM|(c)= - ^|rfx{[6./(4)+a(<. + l)./(<.2)][(</)r(-,c))2-(rf^n}- 

- ^1 d.{[4a(a + l)./(4) + 6a2./(<T2)][(,^r(-,c))2-(4'^)^]} . 

Therefore, the one-loop TK2 kink mass and the semi-classical kink energy are the divergent quantities: 
AMx(c) = AM^(c) + AM|(c) , Esi^^'^ic) = £{$^"^{0)) + /:\Mk{c) 

6 The TK2 kink heat kernel and generaUzed zeta function 
6.1 Zeta function regularization 

We regularize the ultraviolet divergent TK2 kink and vacuum energies in terms of their generalized 
zeta functions: 

AM^(c,.) = ^{^Y^,{C^{c){s)-CkAs)) ■ 

Here, s is the non-dimensional complex parameter already introduced in the previous Lecture; the 
auxiliary parameter of dimensions has also been used in the previous Lecture to keep the 
dimensions in order in the regularization procedure, and the generalized zeta functions are the series: 

N{c) -[^ 1 11 

Therefore, 

AM^(c) = ^lim^ AM^(c,.) = ^ {ck{c){-\) - QkA-\) 

and the divergences reappear at s = — ^, which is a pole of the meromorphic function AM^(c, s) of 
the complex parameter s. 

AM^(c) can also be regularized in terms of generalized zeta functions. Both divergent integrals 
/(4) and /(o"^) when the system is considered in the interval / = — [^^, become the divergent 
series 

Thus, they are obtained as the limits: 



1 ^lL T(s) \m?J ^^0 ^ ' ^ ' u 
2^ _ 1 r(. + l) /^^V^' . , .^ r.22_ , 2 



I(a^) = - lim — • ' • -(^22(5 + 1) , K^^ = -— + a 

^ ' ^lL r(s) y ° f^2;2 
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The regularized contribution of the mass renormahzation counter-terms to the kink mass is: 



AM|(c,s) = hm 



h 



r(s + i) fn 



2 \ 



2L sT-i r(s) 

[6C;,n(s + l)+a(a + l)CK22(s + l)] / dx \ {q^J^^ (x , c)f - {cP^^^ 



+ 



[Ma + l)Cx"(s + l)+6a\,^22(s + l)] / dx {4>^'^ix,c)) 



Therefore, 

AM^(c) = Hm AM|(c,s) 



AL ^[6Ci.ii(^)+^(^ + l)Cx-(^)] j dx [(<^™(x,c))2-(c/.;^))2 
+ [4a(a + l)C^n(i)+6a2c^g.(l)] ( dx [(^^C^, c))^ - 



+ 



6.2 The Cahill-Comtet-Glauber (CCG) exact formula 

The goal in this sub-Section is to compute the one-loop mass shift for the one-component TKl topo- 
logical kink: 



6™^(a;;c = — oo) = -tanhx 



h™{x;c = -oo) = 



The Kink fluctuation operator -K{—(X)) = K for short- is diagonal for any value of a: 



K 



K^^ 




.j4 + 4 6 





dx'^ 





— ^ — 

cosh 2: 



Both diagonal entries K^^ and K^^ are Posch- Teller Schrodinger operators with very well known 
spectra. 

I. The spectrum of K^^: 

There are three types of eigenfunctions 



1. Bound states: 



£0 = 

£1 = 3 



^^0(3;) 

il)^{x) 



1 



cosh X 
sinhx 

cosh'^x 



2. Scattering states: 
with phase shifts: 

3. Half 'hound state: 



i/jkix) = e'''^P2(tanh3;; k) , P2{z; k) = 2,z^ -l-Zikz-k^ 
3k 



6i{k) = — 2arctan 



2-k^ 



£1=4 

2 



tp{k=o){x) = P2(tanhx;0) 
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II. The spectrum of K^'^: 

1. Bound states: I = 0,1,2, ■ ■ ■ , N, N = I[a] 

ei = {2a-l)l , iPi{x) = ^^^2^y_i 2Fi[-l, 2a - l,a - I + + tanhx)] 

2. Scattering states: 

e = q^ + a'^ , ijg{x) = e*^^ 2-F1 [-cr, a,l-iq 



+ e 



Prom these wave functions, one reads the following transmission and reflection scattering coef- 
ficients: 

^^^^ ^ r(cr + 1 - iq)T{-a - iq) ^^^^ _ r(cr + 1 - iq)T{-a - iq)T{iq) 



T{l-iq)T{-iq) ' r(l + a)r(-a)rHg) 

Henceforth, the phase shifts 

identify these scattering processes. 

If (T = E N is a natural number, R{q) = 0, 62 (q) = ^2(1)^ total phase shift 

Hq) ^ larctan f ^-nt-o;(.- - (iV - n)^ + 2.,(iV - .))\ 

3. Half-bound state: If cr = /[cr] = iV G N 

ei = iV^ , ■0^=0(2;) = 2-^1 1; ^(1 + tanhx)] 

2 ^ 

also belongs to the spectrum. 

Thus, one expects that the one-loop TKl mass shift can be calculated exactly from this complete 
spectral information. We distinguish, however, two different situations according to whether or not a 
is a (positive) integer. 

6.2.1 cr = G N: one-loop TKl mass shift from bound states 

If cr is a natural number, the reflection scattering coefficient is zero for both K^^ and K^^. The 
Cahill-Comtet-Glauber (CCG) formula can be applied. This formula gives the one-loop mass shift 
of one-dimensional solitons from the energies of their bound states. Applied to the TKl kink of the 
BNRT model it reads: 



diCosOi) + ^ A^(sina/ — a/cosa/) . (11) 

^ Vj=o /=o / 

11 ^^A r^22. 



The angles are defined in terms of the eigenvalues of the bound states of K and K : 



TT TT ^ ^{2N-l)l ^ 

dn = arccosi- = — , 01 = arccos = — , a; = arccosi-^^ 

u y^i 2 ^26 N 

Thus, for the first five cases we obtain: 
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1. a = 1 



IT 



sin6'o = sinao = 1 



2. a = 2 



3. £7 = 3 



AM{(f'^'^^) = (-- + -^)hm - -hm = -0.984564;im 



1 



1 



vr 2^/3 



AM{cf^^^) = {-- + -^)hm - ( 



1 



vr 2^/3 



ai = arccos( 



smai 



vr 2^3 
V2 



)hni = — 1.33251/i?n 



02 = arccos(2 



sina2 



AM((^^^M = (--+^)nm-(--i(v/5arccos(— )+2v/2arccos(2— = -1.75076/im 
vr 2v3 vr TT 3 3 



4. <T = 4 



ai = arccosi 



smai = - 



3 
4 

15. 



as = arccos( 



.^V3 . 2 

02 = arccos(2^^j , sma2 = - 

1 

sinoa = - 



AM{$^^^) = (-- + ^)^m- (— -i(y7arccos(^) + 2y3arccos(^) + 



,2^3, 



vr 2^/3 



vr vr 



+ ^arccos(-^ )))/im = -2.24628;im 



5. a = 5 



«! = arccosi 



as = arccos( 



V9, 



smai 



smas = — 
5 



a2 = arccos(2 — j , sma2 = - 
5 5 



a4 = arccosi 



,2V6, 



1 

sina4 = - 
5 



TKU 



.3 1 ,^ 

( 1 ^)hm 

vr 2^/3 



15 1, /- ,\/9, ^ ^ ,2^/4, 

( v9arccos( j + 2v4arccos( j + 

TT TT 5 5 



+ \/2Tarccos( ) + 2\/6arccos( )))hm 

5 5 



-2.82180?im 



6.2.2 cT G IR+: one-loop TKl mass shift from zeta functions 

The partition and generahzed zeta functions for the vacuum operator Kq are at the i? — > oo hmit 
respectively!! 



Tre" 



Cxo(s) 



mL 
'2^ 
mL 



dke-''^'+^^+l dqe 



oo 

1 



mL 



4^ r(s) 



mL 



22S-1 



+ 



1 



2s-l 



a 



r(.) 



^In Appendix I it has been shown how this limit can be safely taken, leaving no remnants, when PBC are chosen 
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The poles of Qkq{^) are thus the poles of the Euler Gamma function r(s— i): s— | = 0, —1, —2, • • • , — n, • • •. 
The vacuum energy reads: 



mL 



1 



22S-1 



+ 



1 



r2s-l 



ns-\) 
ns) 



The partition function for the kink operator K accounts for the half-bound state in a subtle way: 



Tr*e-^^ = Tr*e-^^" + Tr^e"^^'' = Tre-'^^'" + e-^^ + - 



+ 00 







N~l 

E 



g-;(2cT-/)/3 _^ ^g-Af(2<7-Ar)/3 _^ 1 
2 vr 







dg— - — e 



and the weight of the last bound state is H = 1, if it is buried in the threshold of the continuous 
spectrum; I[(t] = = cr, or H = 0, if it is not: I[a\ / a. 

Accordingly, via the Mellin transform, the generalized zeta function of the kink operator K reads: 



Ck{s) = CKois)+-+ ' 



N-l 



1 



+- 



-+ 



dk 7r(A;2 + 4)^ /s(2fj - 2iV^(2cj - iV)^ 7o 



d52{q) 



dq 7r(g2 + a'^Y 
(12) 

Both in the partition function and in the generalized zeta function half-zero bound states contribute 
half as much as bound states. 

Because in the general case when a ^ N analytical expressions for the TKl kink generalized zeta 
function are not available, we shall rely on the DHN procedure, very well tested in the paradigmatic 
\<j)^ and sine-Gordon kinks. In this framework, the TKl kink Casimir energy, where the contribution 
of the half-bound state of the vacuum operator Kq - ^ - is always present, can be written as: 



TKlN 



AE{(I) 



+ 



hm 

~Y 

hm 



Ck{s) - Cko{s) - 
Ik^-^V^^ 



TT 



2(2 + 1) 
2(2 + 1 

TT 



2s 



+ 



TT 



(T{a + l) 



2s 



v/3+ir 

Jo 

'N~l 

VW^) + ^ VN{2a - AT) - I + 1 f 
. 1=0 •'^ 



1 



dq (g2 + 



vr 



is still divergent. Zero-point vacuum energy renormalization is not enough. Note that we have sub- 
tracted a finite piece to use the mode-number regularization method. 

The contribution to the one- loop TKl kink mass, induced by the mass renormalization counter- 
terms through the Lagrangian density Cq.t.^ is: 



AM§{^^^) = -2hm [6/(4) + a{a + l)I{a^] j dx (^{(pj^^fi 
= hm [6/(4) + a{a + l)Iicr^)] = hm 



Jo (fc2+4)t 



a{a + l) ^ 1 



27r 



The divergent integrals in AM^{(p^^) can also be regularized by means of zeta function methods: 



/(4) = lim 



2//L 



m^ 



lim • 



m^ 



'^^{k^ + ^Y 
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/(fj^) = lim -i- • • Ci^22(s) = lim • / dq ^ 

We have used the second choice of zeta function regularization proposed in Section §4 because the 
difference in finite renormahzation is included here in AM^. Therefore, the regularized induced energy 
is: 



hm 








hm 


(}^ 







2fiL 

_ 1 

2 



dk 



+ 



dq 



(q^+a^) 



and AM|(<^^i;i 



In the next Table exact results are shown for several values of a, obtained through numerical 
integration of the above formulas. 



a 


AM(TKl)/nm 


0.4 


-0.799335 


0.5 


-0.829892 


0.6 


-0.860369 


0.7 


-0.890955 


0.8 


-0.921788 


0.9 


-0.952966 


0.99 


-0.981384 


1.00 


-0.984565 


1.01 


-0.98775 


1.1 


-1.01664 


1.2 


-1.04925 


1.3 


-1.08242 



cr 


AM(TKl)/fim 




a 


AM(TKl)/;im 


1.4 


-1.11618 




2.4 


-1.4907 


1.5 


-1.15057 




2.5 


-1.53212 


1.6 


-1.18559 




2.6 


-1.57427 


1.7 


-1.22128 




2.7 


-1.61717 


1.8 


-1.25765 




2.8 


-1.65316 


1.9 


-1.2947 




2.9 


-1.70527 


1.99 


-1.32865 




2.99 


-1.74592 


2.0 


-1.33251 




3.0 


-1.75077 


2.01 


-1.33627 




3.01 


-1.75503 


2.1 


-1.37094 




3.1 


-1.79644 


2.2 


-1.41013 




3.2 


-1.84319 


2.3 


-1.45005 




3.3 


-1.89071 



Departure from the reflectionless case (captured by the CCG formula) is thus measured. 
6.3 The high-temperature expansion of the partition function 

For any other TK2 kink, knowledge of the spectrum of the kink fluctuation operator K(c) is grossly 
insufficient to compute the generalized zeta function. Thus, we need to use asymptotic methods to 
obtain sufficiently good approximations to kink generalized zeta functions. 
The heat equation kernel for the i^Q-heat equation 



a 



80 doP 





+ 4 



a 





a^ 



a/3 



KKo{x,y;(3)=0 , KKo{x,y;0) 



6{x - y) 






5{x - y) 



of the vacuum fluctuation operator Kq^ for small (3 is: 



dx' 



+ 4 











KKo{x,y;P) 



4/3 (=^-v) 

e 



4/3 



V 



g 47r/3 



J 



The TK2 kink fluctuation operator is the 2x2 matrix Schrodinger differential operator 



K{c) 



■i& + 4 + y"(x;c) 



Ko + V{x-c) 



dx' 
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whereas the corresponding heat equation kernel 



^ (9/3 + ^ / J 



with 



6{x - y) 
6{x - y) 



can be written in the form 

KK{c){x,y;P) = CK{x,y;P) ■ KKo{x,y;f3) , Cx(a;,y;0) 
if the 2x2 matrix CK{x,y; P) satisfies the transfer equation 



1 
1 



d X — y d d"^ 
dp ^ p dx 



and if it is set to the unit matrix at infinite temperature. 
Solving the transfer equation as a power series in P 

oo 
n=0 



4''ix,y;K)=6 



AB 



the PDE transfer system of equations becomes tantamount to the system of recurrence relations: 



in+l)c^^,{x,y;K) + {x - y)-^{x,y; K) 



- (4(^^^ -6^')+ a\5^' - 5^'))^ (X, y; K) - g V^^^(x)c^^(x, y- K) 



The matrix elements of the heat equation kernel for K{c) have the high-temperature asymptotic 
form 



p-4/3 



e 4/3 



2 OO 



(x~y) 



2 OO 



n=0 



K%){x,y^P) 



4/3 



■e-^^-Y,c^n{x,y;K)P^ 



n=0 



•e 4/3 



n=0 



n=0 



Actually, we are interested in the trace of the heat kernel, both matrix and functional, to find the 
partition function for small p. The recurrence relations become 



1 



n+l 



C=l 



when y — > X. To deal with this delicate point, we have introduced the following notation: 

^AB 



«C„^^(x) = linr^(.,,;K) 



f^k AB 

^^^C^^{x) = lim '^-^{x,y-K) = 6^'5^^ 
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We also need (obtained after differentiating the first recurrence formula /c-times) recurrence relations 
among the derivatives: 



The high temperature asymptotic expansion of the partition function reads: 



^/47^p ^ — ' L J L— +00 / mL 



n=0 



2 



L — ^OO / mL 

2 



Using the recurrence relations, the c^^{x, x; K) densities can be found. They are the conserved charges 
of a generalized matrix KdV equation, see Appendix III. 



6.4 The Mellin transform of the asymptotic expansion 



We now write the generalized zeta function of Kq as the Mellin transform of the Kq partition function 
split into two integrals: 



mL 1 



\/4^ r(s) 

mL 1 



+ 



^ r(s) \4s-\ 



7[s-i 4]+r[s-i 4] 



dj3l5' 
1 



+ 



r2s-l 



The incomplete ^[s — ^, 4] and ^[s — ^, o"^] have poles at s — ^ = 0, —1, —2, —3, • • • but their comple- 
mentary functions V[s — ^,4] and T[s — ^,cr^] are entire functions of s. 

To obtain the generalized zeta function from the asymptotic expansion of the K partition function 
the Mellin transform is also split into two integrals, inside and outside the convergence radius: 



cm 
1 



T{s) 

oo 

V4^ „J77^, 1 \ 



n=0 



4' 



s+n- 



n=7V()+l 



s+n- 



(T 



2s+2n-l 



+ 



1 



47r 



^2s+2n-l 



+ 



The two zero modes have been not accounted for and the incomplete Gamma functions ^[s + n— ^,4] 
and ^[s + n — ^,ct^] have poles at s-\-n— ^ = 0, —1, —2, —3, • • •. A large but finite number Nq is chosen 
to separate the contribution of the high-order coefficients 



n=Afo+l 



AS+n- 



which are holomorphic functions of s for Res > —Nq — \. Bk{s) = df3 P'"^ ^Tr*e , however. 



is a entire function of s. 
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6.5 The high-temperature one-loop TK2 kink mass shift formula 

Neglecting the (very small) contribution of the entire functions, the TK2 kink Casimir energy becomes: 



AM 



c 



. 1 V 

2 



1 



No 



7[s + n- i 4] 



+ cf{K{c)) 



n=l 



7[g + n- 

fj2s+2n-l 



i.e., the zero-point vacuum energy renormalization takes care of the term coming from Cq^{K{c)) and 

The other correction due to the mass renormalization counter-terms can also be arranged into 
meromorphic and entire parts: 



2V47r 



1 



' r(.) 



+ 



4^+5 

c?{K{c)) 



7[s + ^,4]+r[s + ^,4] 



+ 



The mass renormalization terms exactly cancel the contributions of c\^{K{c)) and cf^{K{c)). Our 
minimal subtraction scheme fits in with the following renormalization prescription: in theories with 
only massive fluctuations quantum corrections vanish at the limit where all the masses go to infinity. 
We end with the high-temperature one-loop TK2 kink mass shift formula: 



AMx(c) 



hm 



In this case, the subtraction of the two zero modes contributes to the mass shift in the c-independent 
quantity: 

(0) _ hm 



AM 



K(c) 



-0.56419nm 



TT 



i.e., for each kink in the TK2 family we must subtract the same quantity to discard zero mode effects 
at the one-loop level. 



6.6 Mathematica calculations 

Computational limitations put a practical bound on the choice of Nq. Knowledge of, say, ^^^C2 requires 
computation of 36 (or 9A^^ in field theories of N scalar fields) densities: 

WCo (3) Co (2) Co (i)Co WCo 
(2)Ci 

(0)C2 

In general, evaluation of ^'^^Cnix) requires previous calculation of 

4(1+3+5+7+- • •+2n-l+2n+l) = 4(n+l)2 , Ar2(i+3+5+7+. . .+2?i-l+2n+l) = N'^{n+lf 

This count could be slightly abbreviated bearing in mind that the 4(2?i + 1) coefficients in the upper 
row are fixed by the initial conditions of the recurrence relation. 
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6.6.1 TKl topological kinks 



The previous formulas can be applied to the c = — oo case for several values of u; i.e., to compute, by 
means of the asymptotic method, the TKl kink mass shift, to find -using Mathematica- the results 
shown in the next Table. 



a 


AMK(-oo)/^m 




a 




0.5 


-0.962386 




1.1 


-1.05073 


0.6 


-0.970537 




1.2 


-1.07468 


0.7 


-0.981183 




1.3 


-1.10097 


0.8 


-0.994487 




1.4 


-1.12939 


0.9 


-1.01053 




1.5 


-1.15971 


1.0 


-1.0293 




1.6 


-1.19174 



a 




1.7 


-1.22526 


1.8 


-1.2599 


1.9 


-1.29571 


2.0 


-1.33324 


2.1 


-1.37074 


2.2 


-1.41007 



By comparing these numbers with those obtained by the exact procedures of Cahill-Comtet- 
Glauber and Dashen-Hasslacher-Neveu we show that the error of the asymptotic method decrease 
with increasing a, see next Table and Figure: 



Value of a 


a = 0.9 


a = 1.5 


a = 2.0 


a = 2.2 


Relative Error 


6.0 % 


0.79 % 


0.055 % 


0.004 % 



-0.4r 

-0.6 
-0.8 



-1.2 
-1.4 
-1.6 
-1.8 



AM 



K 



O 

gpnn'% . 

. 5 5^0.1 . 5 2 2T5 3 375 



Figure 6: One-loop correction to the one-component topological kink (TKl) mass in units of Km. •, DHN 
formula . asymptotic series. 



6.6.2 Elliptic two-component TK2 topological kinks 



Along as cr < 1, c = is a kink orbit. This kink orbit is a half-ellipse in the (f)i- 
corresponding kinks are TK2 kinks for which the kink fluctuation operator reads: 



b2 plane and the 



Km 



+ 4- 



2(l-f2a^) 
cosh^[2(l-a-)a; 
sinh]2(l-o-)x 



^2 



sinh[2(l-o-)a] 
cosh^[2(l-CT)a;] 

cr(5-7o-) 
cosh^[2{l-(7)x] 



The ri^2(0) Seeley coefficients are computed for several values of a in the a € [0.96, 1) range. Together 
with the one- loop TK2(0) mass quantum corrections according to the asymptotic method, they are 
shown in the next Tables: 







0.96 




0.97 


n 






C(^(0)) 


eft {Km 


1 


12.6806 


3.83333 


12.5025 


3.87629 


2 


27.2626 


2.91223 


26.3923 


2.84445 


3 


43.0467 


0.929709 


40.8995 


0.971713 


4 


51.8330 


0.431151 


48.3164 


0.391887 


5 


49.4842 


-0.00452478 


45.2317 


0.0193729 


6 


38.8786 


0.0514525 


34.8351 


0.0380063 


7 


25.8993 


-0.0156285 


22.7417 


-0.0090224 


8 


14.9659 


0.00713977 


12.8764 


0.00456185 
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0.98 


a = 


= 0.99 


71 






V-'*- J 


(K(Q)) 

V-**- J 


1 


12.3299 


3.91837 


12.1624 


3.9596 


2 


25.5597 


2.78109 


24.7630 


2.7219 


3 


38.8820 


1.00819 


36.9849 


1.03968 


4 


45.0733 


0.358188 


42.0798 


0.329352 


5 


41.3849 


0.0389505 


37.9012 


0.0548776 


6 


31.2486 


0.0273206 


28.0632 


0.0188943 


7 


19.9962 


-0.00450157 


17.6055 


-0.000842262 


8 


11.0960 


0.00265278 


9.57626 


0.0003676 



a 




0.96 


-1.06082 


0.97 


-1.05253 


0.98 


-1.04422 


0.99 


-1.03624 



6.6.3 One-loop breaking of classical kink degeneracy 

To end this Section we compute one-loop mass shifts for TK2 kink families by means of the asymptotic 
approximation. In the next Tables we offer results for a = 1.5, a = 2, and a = 2.5 and several values 
of c between —30 and a value very close to . 



a = 1.5 




cr = 2.0 




c 


AA/ 


c 


AM 


-30 


-1.16009 


-30 


-1.33281 


-27.5 


-1.16017 


-27.5 


-1.33281 


-25 


-1.16128 


-25 


-1.33281 


-22.5 


-1.16042 


-22.5 


-1.33281 


-20 


-1.16061 


-20 


-1.33281 


-17.5 


-1.16088 


-17.5 


-1.33281 


-15 


-1.16128 


-15 


-1.33281 


-12.5 


-1.16193 


-12.5 


-1.33281 


-10 


-1.16313 


-10 


-1.33281 


-7.5 


-1.16597 


-7.5 


-1.33281 


-5 


-1.18205 


-5 


-1.33280 


-4.6801886 


-1.24345 


-4.001 


-1.33280 


-4.68018860186678332 


-1.25103 


-4.00001 


-1.33280 



cr = 2.5 


c 


Ail/ 


-30 


-1.52784 


-27.5 


-1.52782 


-25 


-1.52780 


-22.5 


-1.52778 


-20 


-1.52774 


-17.5 


-1.52769 


-15 


-1.52760 


-12.5 


-1.52744 


-10 


-1.52711 


-7.5 


-1.52626 


-5 


-1.52285 


-4 


-1.52168 


-3.97 


-1.52915 


-3.96594571 


-1.55402 


-3.96594570565808127 


-1.56127 



-30 



-25 



-20 



-15 



-10 



a=l . 5 



a = 2 . 



-0 . 9 

-1 . 1 
-'l.2^ 

-1 . 3 



AAiAiAiA' 

a = 2.5 



The One-loop Quantum 
Mass Correction in the cases 
a = 1.5, a = 2.0 and a = 2.5 



1 . 4 
1 . 5 



1.6 



Am 
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The behavior is the same for a = 1.5 and a = 2.5: the classical degeneracy of kink energies 
survives one-loop quantum fluctuations for values of c lower than the critical values where TK2 kinks 
start to split into two lumps. The one-loop mass shifts for TK2 kinks formed by two lumps are 
remarkably higher (in absolute value) and increase with kink separation. For the value cr = 2, kink 
energy degeneracy is not lifted by one-loop fluctuations. The reason is that the BNRT model for this 
value of a is no more than two A(/)^ independent models if appropriate linear combinations of (pi and 
02 are chosen. 



7 The planar Abelian Higgs model 

Given a complex scalar field and a U{1) gauge potential 



d 



B^{y^)Tr- ■ — ' LieU(l) 



the action for the planar Abelian Higgs model reads: 



2\2 



d 



G 



where the volume and the metric tensor in (2+l)-dimensional M^'^ Minkowski space are given below, 
together with the dimensions of the fields and parameters: 



d% = dy^dy^dy^ 
Defining non-dimensional coordinates, fields, and parameters. 



Qi^y = diag(l,-l,-l) 
e] = [A^] = Ar^L^^ 



y^ = —x^" , ip = v(j) = v{(pi + i(j)2) , Bf, = vAf, , = A 



ev 



the action and the field equations are: 



S 



d-^x 



\f,.F^' + \{D,<prD^ct> - ^{<P*ct>{xn - If 



D„ 



F, 



dxt" dx" 



Df,Df'(j) = —m 



7.1 Feynman rules in the Feynman-'t Hooft R-gauge 

There is C/(l)-gauge symmetry 

and if A G [0, 2tt] is a constant angle the vacuum orbit of the gauge group is: 



da 



A]:{xn = o, , K 



da 
dxi^ 



(xn 
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We shift the scalar field away from the vacuum in H{x^)-iliggs and G(x'')-Goldstone fields: 

(l>(x'') = l + H{x>') + iG{x^') = 
The choice of the Feynman-'t Hooft R-gauge 



(1 + Hy{x^') = cosa(x^)(l + H{x>')) - sma{x^')G{x>') 
G'ixi") = sina(x'^)(l + H{x^')) + cosq(x'^)G(x'^) 



R{A^,G) = d^A^{x^')-G{x^ 



needs a Faddeev-Popov determinant to restore unitarity, which amounts to introducing a complex 
ghost field. Because 



we find 



Det- 



5R 
5a 



R{A'^, G') ~ RiA^, G) + (□ - 1 - Hixf")) ■ 5a(x^) 
[dx*{x>^)][dxixn]ew{iS^host[x\x]) 

[dx*ixn][dx{xn]^xp{i [ d'xx*ixnia-i-Hixn)x{xn 



All this together allows us to write the action in the form 



S + Sfff, + 5, 



g.f. ~r 'Jghost 



1 



-^A,[-g>^^n + 1)]A, + d,x*d>'x - X*X 



1 



1 



+ -dMd^'G - -G^ + -d^Hd^'H 



k2 



-H{H^ + G^) + A^id'^HG - df^GH) + H{Af,A'' - x*x) 



— {R' + Gy + -(G" + R')A^,A>' 

o 2, 

which encodes the Feynman rules shown in Tables 7 and 8. It should be noted that fermionic (ghost) 
loops carry a (—1) factor. 



7.2 Plane waves and vacuum energy 
7.2.1 Vector bosons 

The general solution of the linearized equation for small fluctuations of the vector field around the 
vacuum 

^ - VV + Ij M^(xo,x) = , A^(xo,x) ~ 0^ + M^(xo,x) 



is the plane wave expansion 



bA^{xQ,x) 



EE 



2uj{k) 



if the dispersion relation /cq — /c/c — 1 = holds. Here, we denote 



kx knX^ /Cq37q JvX 



kk + 1 
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Table 7: Propagators 



Particle 



Field Propagator Diagram 



Higgs 

Goldstone 

Ghost 



Hix) 
G{x) 



ieh 





- + ie) 




ieh 




-1 + ie) 




ieh 


v{k^ 


-1 + ie) 




iehgt"" 



Table 8: Third- and fourth-order vertices 



Vertex Weight 



Vertex Weight 




2 ^ 
-61K — 

he 




-IK 



he 



he 



he 



2 ''^ 
— 

he 



-Sin — 



he 



-IK 



he 




he 



2i^g^^^ 
he 
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and consider periodic boundary conditions on a square of area nn?L?', m = ev, such that: 

ki = — -rii , nj e Z ,z = l,2 
mL 

The polarization vectors eJJ(/E), a = 0, 1, 2, satisfy the ortho-normahty condition: 

e"(A:) • ef^ik) = e"^(A;)e^(yfc) = -(-I)''"" J"'^ 
The hnear classical Hamiltonian 



m] = y / d^x5A^{xQ,x) 



dxl 

\<5qO 



+ VV- 1 



5A^{xo,x) 



{-lf-'^{al{k)ao^{k) + ao,{k)al{k)) 



leads, via canonical quantization 

[a^{k),aim = {-^f'''5ap5^u^ 
to the quantum Hamiltonian for free massive vector bosons: 



The contribution to the vacuum energy of the massive vector bosons is: 



= E E ^'^(^) = ^TV[- VV + 1] 



7.2.2 Higgs bosons 

The linearized field equations for small fluctuations of the Higgs field around the vacuum, 
are solved by the plane wave Higgs expansion 



6H{xo,x) = 



SH{xo,x) = —S2 
vL V ev ^-^ 



2uj{k) 



a*{k)e'^'' + a{k)e 



-ikx 



io{k) = +y kk + Hi 



The classical energy of the Higgs plane waves 



d6H 

dxQ 



+ V6HV6H + K^5H6H 



hm 



^uj{k) a* {k)a(k) + a(k)a* (k) 



becomes -through the canonical quantization [a{k) , a"^ (q)] = 6^^- the quantum Hamiltonian for the 
Higgs bosons: 

H^'^\5H] = /im^cj(fc) (^\k)a{k) + ^ 
The contribution to the vacuum energy of the Higgs bosons is thus: 

= E ^^(^) = ^Tr[-VV + 2 
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7.2.3 Goldstone bosons 



Simili modo, the linearized field equations for small fluctuations of the Goldstone field around the 
vacuum ^ 



are solved in terms of Goldstone plane waves: 

SG{xo,x) = —"S^- 
vL V ev ^-^ 



a*(Ai)e*^^ + a(A:)e- 



-ikx 



k^-kk -1 = 



The classical energy of the Goldstone plane waves 



io{k) = +ykk + 1 



fd5G_ 



+ V5GV5G + 5G5G 



hm 



a* {k)a{k) + a{k)a* (k) 



is promoted through canonical quantization [d{k), a{q)] = (5^^ to the quantum free Hamiltonian 

H^^^[6H] = hmY^uj{k) (a^k)a{k) + ^ 



such that the contribution to the vacuum energy of the Goldstone bosons is: 

= E = ^Tr[-VV + 1] 2 . 



7.2.4 Ghost particles 

The contribution of ghosts to the vacuum energy is more tricky. The solution of the linearized field 
equations for small fluctuations of the Ghost fleld around the vacuum 



dx, 



2 VV + 1 ) 5x(xo,x) =0 



is also a plane wave (ghost) expansion: 



1 / h 



vL V ev ^-^ 



2uj{k) 



c(A;)e~^^^ + d*(A;)e 



ikx 



k^-kk -1 = 



uj{k) =+\lkk + l 



The coefficients, however, are Grassman variables -classical cousins of Fermi fields- satisfying the 
anti-commutation relations: 



c^{k) = (f{k)=0,yk , c{k)d{q) + d{q)c{k) = 

c{k)c{q) + c{q)c{k) = d{k)d{q) + d{k)d{k) = , Vfc, g 
c{k)c* (q) + c* {q)c{k) = d{k)d* (q) + d* {k)d{k) = 
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The classical energy of ghost plane waves looks familiar up to a sign 

'dSx* d5x 



dxQ dxQ 



= ^^uj(k) \c* {k)c{k) + d* (k)d{k) - c{k)c* (k) - d{k)d* (k) 
k 

but canonical quantization proceeds by the anti-commutators 

{cHk),c{q)} = {dHk)Jm=%^ , 
and the free quantum Hamiltonian is: 

H^^'^[6x] = hm'^ujik) (^ct(fe)c(fe) + d}(k)d(k) - 1 

k 

Thus, the contribution to the vacuum energy of Ghosts is negative: 

AE^i^^ = - ^ hmuj{k) = -n?nTr[-VV + 1]5 

k 

Preserving unitarity, the ghosts exactly cancel the contribution to the vacuum energy of the non 
physical temporal vector bosons and Goldstone bosons. 
In sum, the vacuum energy in the planar AHM 



r=l 



/^Eq = ^ A^J''^ = /imTr[-VV + 1]5 + ^Tr[-VV + k^]^ 



is due only to Higgs particles and transverse massive vector bosons, as it should be. 

7.3 One-loop mass renormalization counter-terms 

Denoting as I(c^) the divergent integral 

d'^k i 



(27r)3 fc2 - c2 + ie 
the AHM encompasses the following one-loop divergent graphs: 
1. Higgs boson tadpole: 



+ \ ) + ' ; + — 

= -2i{K^ + 1)/(1) + finite part 




2. Higgs boson self-energy: 




-2z(k^ + 1)/(1) + finite part 
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3. Goldstone boson self-energy: 




= -2i{K^ + 1)/(1) + finite part 



4. Vector boson self-energy: 

-|- [ ] -\- ■\/\f\C %vAA/ 

= 2il{l)g'"' + finite part 

These calculations are performed in a detailed manner in Appendix §. 12.4. 

All the finite parts are proportional to 1{k?) — 1(1) so that with a minimal substraction scheme we 
can get rid off all these one-loop divergences by adding the counter-terms 

rf,. = h{K' + 1)1(1) [\^\' - 1] , Lit, = -him^A'^ 

to the Lagrangian. 

In (2-|-l)-dimensions however, the graphs above are the only divergent diagrams in the system 
(the theory is super-renormalizable). Thus, the diagrams coming from these counter-terms completely 
cancel any divergence (not only to one-loop order) arising in the vacuum sector of the model: 




Diagram Weight 



X 


2i{ 


y + 1)^(1) 


X — 


2i{ 


y + 


-------------X — 


2i{ 


y + 1)^1) 






-2il{l)g^^ 



We stress, however, that renormalization in the planar Abelian Higgs model requires more than 
merely normal ordering, contrarily to (l-l-l)-dimensional scalar field theory, where normal order is 
sufficient to cope with ultraviolet divergences. In this system there are divergences due to graphs with 
two vertices that are not suppressed by normal ordering. 

7.4 Self-dual Abrikosov-Nielsen-Olesen vortices 

There are also topological solitons in the AHM: solutions of the static field equations 

d^F,j = Jj ■ DiD,(t>= —(I}{4>*(t>-1) 

where Jj = | {(j)*Dj(j) — {Dj(j))* (p) is the electric current, of finite energy: 

E = I d^x[^FijF,, + ^{D,^)*Di<p + '^{<p*<p - 1)2] . 
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In fact, the configuration space C 

C = {(j){x) e Maps{R'^X),Mx) G Maps{R^,TR^)/E{<l),Ai) < +00} 

is topologically disconnected. The boundary of the spatial plane is the sphere of infinite radius 
= limr_>_|_oo [xf + X2 = r^] . Finite energy configurations comply with the asymptotic behavior 

4>*4>\si^ = ^ , D,cp\si^ = id,(t)-iAi4>)\si^=0 , 

such that 

9 = arctan^ , cP\si^ = e'^^ ,l£Z , Ai\si^ = -icj)*di(P\si^ 

provides the map lim^^+oo (t'i^) '■ S^^ — > Sl between the sphere at infinity and the vacuum orbit §1. 
Continuous maps between one-dimensional spheres are classified according to the first homotopy group 
and, because the temporal evolution is continuous, Ilo{C) = Ili{Sl) = Z, the zero homotopy group 
of C is non-trivial. Thus, C = U^g^C; is the union of disconnected sectors characterized by an integer 
number I and the magnetic flux of any finite energy configuration is quantized: g = J d'^xFi2 = 27rL 

We shall restrict ourselves to the critical point between Type I and Type II superconductivity: 
= 1. The energy can be arranged in a Bogomolny splitting: 

E = 1^ {\Di<P ± iD^cPl' + [F12 ± i(</.> - + \\g\ . 

We immediately realize that the solutions of the first-order equations 

Di<P±iD2(l) = ; Fi2±^((/.*</.-l) =0 

are absolute minima of the energy, and are hence stable, in each topological sector with a classical 
mass proportional to the magnetic flux. 

Because the vector held is asymptotically purely vorticial, these solitonic solutions were christened 
as vortices by their discoverers Abrikosov, Nielsen and Olesen. Also, since the first-order equations 
can be derived from the self-duality equations of Euclidean 4D gauge theory through dimensional 
reduction, the ANO vortices are called self-dual at the limit = 1. 

7.5 Self-dual vortices with spherical symmetry 

Another simplification is to consider the spherically symmetric ansatz: 

(j)i{xi,X2) = fir)cosW , (/)2(a:;i,X2) = /(r)sinZ6' 

ct(r) ct(r) 

Ai(xi,X2) = —I sinO , A2{xi,X2) = I cosO 

r r 

The first-order equations reduce to 

to be solved together with the boundary conditions 

lim /(r) = 1 , lim a{r) = 1 

r— >oo r— >oo 

/(O) = , a(0) = 

dxiAi = —I (h 5 = 27rZ , 



r=oo 



required by energy finiteness plus regularity at the origin (center of the vortex). 

To solve this non-linear ODE system, we follow the three-step De Vega-Schaposnik procedure 
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1. For small values of r, a power series is tested in the first-order differential equations. 

2. The first-order ODE system is solved exactly for large r: the asymptotic behavior of the solutions 
is thus found. 

3. Finally, a numerical scheme can be implemented by setting a boundary condition at a non- 
singular point of the ODE system, which is obtained from the power series for small values of 
r. The behavior of the vortex solutions for intermediate distances is then described by means of 
an interpolating polynomial. 

In terms of the field profiles f{r), a(r) the magnetic field and the energy density are: 

Bir) = l^ , .(r) = i(l-/W + ^(l-«Mmr) • 

The field profiles, the magnetic field and the energy density are shown in Figure 1 for / = 1,2,3,4. 
B{r) always has a maximum at r = but the origin is only a maximum of the energy density for 
I = 1. 




Figure 1. Plots of the field profiles a{r) (a) and f(r) (b), the magnetic field B{r) (c), and the energy density 
e{r) for self-dual vortices with I = 1 (solid line), I = 2 (broken line), I = 3 (broken- dotted line) and I = 4 

(dotted line). 

In Figure 2, three-dimensional plots of the energy density are shown for / = 1, 2, 3, 4. 




Figure 2. 3D graphics of the energy density for l~l, l^2,l~3 and I ~ 4 self-dual symmetric ANO vortices. 
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Figure 3. (left) The I = l-vortex scalar iso-vector field s(yXi,X2) = Res(a;i,X2)i + lTa.s{xi,X2)j projected on the 

R spatial plane (hedgehog) . 
(right) The I =^ l-vortex vector field V{xi,X2) — Vi{xi, X2)i + V2{xi, X2)j ■ Note that V • V{x) = 
and V A V{x) — k. Also, V{0, 0) — l/(cx)C0s6', cx)sin0) = (Poincare-Hopf theorem) . 



7.6 Two-vortex solutions with different centers 

I = 2 ANO self-dual solutions formed by two / = 1 vortices with centers separated by a distance d can 
also be obtained approximately. We shall implement the variational method of Jacobs and Rebbi in 
two stages: 

1. First, trial functions, depending on a single variational parameter w, 

4>M,z*) = ^z,z*)[u;f^'H\z-d/2\)fW(\z + d/2\) + 

\z^ - {d/2f 



A^{z,z* 



(1) I 



z* - d/2 
2i 



Q^^^(|z-d/2|) + 



z* + d/2 



a«(|z + d/2|) ) + 



+(l-u:)-a^'K\z\ 



z* 

are built. Here z = xi + ix2, A'^{z, z*) = Af{z, z*) + iA'^{z, z*) and 

Hz,z*) 



z^ - {d/2Y 



z*2 - {d/2Y 



g = An 



f^^\ a^^\ /(^^ and a^^^ stand for the functions / and a associated with self-dual solutions with 
cylindrical symmetry, respectively, with vorticity / = 1 and 1 = 2. 

Plugging this ansatz into the energy functional, expression E[lo) is set to be minimized as a 
function of lo. 

2. Then, a deformation is added such that: 1) the scalar field vanishes at the two centers; 2) 
the gauge-invariant quantities associated with the solution are symmetric with respect to the 
reflection z — > z*. 

The invariant ansatz 



{z,z*) = <j)^{z,z*) + 



N i 



+ ^{z, z*) \z^ - {d/2f\ (cosh /^^-^ 



1=0 j=0 



Z \3 I Z 
Z* J \ Z 



* \ J 
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A{z,z*)=A'{z,z*) + 

N i 



+ 



1 



cosh \z\ 



contains H = 3 



i=Q j=0 



{N+l)(N+2) 



Z \1 I Z 

z* ) \ z 



*\J 



N i 



i=0 j=0 



[zz*y 



Z \0 I z 

z* ) \ z 



* \ 1 



variational parameters: fij, af^-, afj. 



A three dimensional plot of two- vortex solutions is shown in the next Figure for distances d = 1,2,3. 




Figure 4. 3D graphics of the energy density for I = 2 self-dual separate vortices with centers at distances d = 1, 
d = 2, (i = 3. Here we have plotted the cruder approximation: N ^ 1. 



7.7 The vortex Casimir energy 
7.7.1 Bosonic fields 

Let us consider small fluctuations around vortices 
(l){xo,x) = s{x) + 6s{xq,x) , 



Ak{xo, x) = Vfc(x) + (5afc(xo, x) 



where by s{x) and Vk{x) we respectively denote the scalar and vector field of the vortex solutions. 
Working in the Weyl/background gauge 

^o(2;o,^) = , dk6ak{xo,x) + S2{x)6si{xo,x) - si{x)5s2{xo,x) = , 

the classical energy up to 0(5^) order is: 



^ + ■'^g.f. + ^ghost 



d x 



^ {xo,x)^ — {xo,x) } + 



\ dxo ^ ^ dxQ 
+y / d^x{S^''{xo,x)K6^{xo,x) + 5x*{x){-A+\s{x)\^)6x{x)} 
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where 



and 



K 



( - A+|s|2 




/ 5ai{xQ,x) \ 
5a-2{xQ,x) 
5si{xQ,x) 
\ 5s2{xq,x) j 





A+|s|2 



A = VV , 



-2ViS2 
-2V2S2 



-2ViS2 -2V2S2 - A+i(3|s|2 + 2yfcyfc-l) 
V 2ViSi 2V2S1 Wkd^, 

The hnearized field equations read 



dxl 



(xo, x) + Kbi{xQ, x) = 



2ViSi \ 
2V2S1 

A+i(3|s|2 + 2yfcVfc-l) j 



A+|s(x)P)5x(^') =0 



or, in more detailed form, 

(- A +|s(f)|^) ba\{x^,x) = IV \S2{x)bs\(xQ,x) - 2\I \S\{x)bs2{xQ,x) 
(- A +|s(f)p) 502(3:0, = 2V2S2(^)5si(xo, - 2\l 2S\{x)bS2{xQ,X) 
A +i(3|s(x)|2 + 214(f) Vfc(f) - 1)) 5si(xo, f) - 214(x)afc<5s2(xo,f) = 

= 2ViS2(a^)5ai(a;o,3f) + 2V2S2(f)(5a2(a;o, x) 
A +^(3|s(f)|2 + 214(f) yfc(f) - 1)) 5s2(xo,f) -2I4(x)afc<5si(xo,f) = 

= — 2ViSi(f)5ai(a:;o, f) — 2V2Si(f)(5a2(a;o, f) 

Let us assume the following spectral resolution of K in the orthogonal complement to its kernel 
(no bound states): 



J]] E:ABuiJ^(f;fc) = e(fc)u^^^(f;A?) ,1 = 1,2,3,4, ei$) = kk^\ 



B = l 

u^l\x]k) = e'''^v^^\x; k) 
4 



lim v^J^{x]k) = VA{k)6^j^ 



[ d'xuf*ix;k)u^^\x;q)=e'v'L\/'' 

A 1 J 



A=\ 

Expanding the small fluctuations in terms of the positive eigenfunctions of K 

1 [fl V - ^ 

vL V ef 



<5Ca(3;o,^) 



1 



a*i(k)e'^^^u^^^\x; k) + a/(A:)e-*"^«n^^(f; k) 



in. 



k 1=1 ^2e{k) 

one ends with the classical free Hamiltonian: 

4 

i/(2) + //(2) =h'^Y.Y. ^(^) [o*i{k)ai{k) + ai{k)a*j{k) 
h 1=1 
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7.7.2 Ghost field 

Assuming also that there are no bound states in the positive spectrum of i^*^ = — A 
K'^u{x;k) = e{k)u{x;k) , e{k) = kk + 1 

u{x; k) = e^^^v{x; k) , lim v{x; k) = v{k) 



Sx u* {x; k)u{x; q) = e^v^LF'd 



k q 



small fluctuations of the ghost field can be expanded in terms of positive eigenfunctions of K*^ 



1 



Sx'a{xo,x) = —r\— ■ 

vL V ev ^-^ 



1 



c{k)u*{x; k) + d*{k)u{x; k) 



k \/Hk) 

with Grassman variables as coefficients. The ghost classical free energy is thus 



(2) 

Ghost 



Yl \c*{k)c{k) + d*{k)d{k) - c{k)c*{k) - d{k)d*{k) 



Two remarks are in order: (1) Because the c's are Grassman coefficients there are no cc or c*c* terms, 
and c*{k)d{—k) + d* {k)c{—k) = 0. (2) Note also that the ghost fields are static in this combined 
Weyl-background gauge. Therefore, their energy is one-half with respect to the time-dependent case. 
The canonical quantization 



[ai{k),d\{q)]=6ij6i:^ 
leads to the quantum free Hamiltonian 

H^'^+H^l+HiL = nm.Y^eCk) 



{c{k),c\q)}=6 



k q 



{d{k),S{q)} = 6^^^ 



^ (d\{k)dj{k) + U\k)c{k) + d\k)d{k) - 1 



.1=1 



such that the vortex Casimir energy reads 



^Ey = ^STr* kI = ^Tr* - ^Tr* {K^f^ 
^ 2 2 2 ^ ^ 



In a similar manner we write the vacuum energy: 



Ai.o = ^STr4 = ^TY4-^Tr(i.o-)^ 



/ -A + 1 



V 





-A + 1 








-A + 1 










A+1 



+ 1 / 



The zero-point vacuum energy renormalization, defining the vortex Casimir energy, is performed in 
the formula 



/\M^ = AEv - AEo 



hm 



STr*K2 - STrK^ 



63 



8 The vortex heat kernel and generahzed zeta function 
8.1 Deformation of the first-order equations 

The dimension of the kernel of the operator K ruhng the small vortex fluctuations orthogonal to the 
gauge group is the dimension of the moduli space of vortex solutions. Small deformations around 
vortices 

4>{x) = s{x) + 5s{x) , = ^ji^) + Saj{x) 

are still solutions of the first-order equations if 



'12 



T^P) —d26ai + di5a2 + siSsi + S26S2 = 
(^202 - A2(j)i) = si6ai - S26a2 - {82 - Vi)Ssi - {di + V2)5s2 = 



(^2(^1+ ^202) + {di(l)2-Ai(pi) = 0<^S26ai + si5a2 + {di + V2)6si-{d2-Vi)6s2 = 

To avoid pure gauge deformations, we set the background gauge: 

dj5aj{x) + S2{x)5si{x) — si{x)5s2{x) = 

Therefore, the tangent space to the moduli space of self-dual vortices is the Kernel of the first-order 
deformation operator P: 





-d2 


di 


Sl 


S2 \ 


1 ^O'l 


[x) \ 




-di 


-82 


-S2 


Sl 


6a2 


ix) 




Sl 


-S2 


-82 + Vi 


-di - V2 


5si 


[x) 


v 


S2 


Sl 


di + V2 


-d2 + Vi 1 


\ 5S2 


[x) ) 



One easily checks that K = V'^T) has a supersymmetric partner: K = VD^ . 



Vi{x) = 



K- 



( -A+|s|2 



V 





A+|s|2 
- A+i(|s|2 + l) + 14yfc 












1YA 



A+i(|s|2 + i) + yfcyfc / 



The index of the deformation operator - ind T) = dimKerP — dimKerP^^ - is in this case equal to the 
dimension of KerX because dimKerP^ = 0, K~ being definite positive. 

8.2 The kernel of the heat equation 

The heat equation kernel of a x matrix differential operator of the general form 

K = KQ + Qk{x)dk + V{x) 
is the solution of the i^-heat equation 

l^^l + K^KK{x,mf3) = Q , 

with initial condition: Kk{x, y;0) = I • 6^'^^ {x — y). From the kernel, one derives the partition function 

Tr e"^^ = tr / d^xKxix, x; (3) 
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which, in turn, provides the dimension of the vortex moduli space 

ind P = Tr e"'^^ - Tr e"^^" 
because Non-zero Specif = Specif". To find the kernel one writes 

Kxix, y; /?) = Ck{x, y; (5)Kkq{x, y; /3) 

where 

KK,{x,mf3) = — -l-e 
is the -RTo-heat equation kernel for small (3. CK{x,y; (3) satisfies the transfer equation 

+ ^^^{dkl - ^Qk) - AI + QA + y| Ck{x, y; /?) = , 

and is the unit matrix Ck{x, y; 0) = I at infinite temperature. 

8.3 The high-temperature expansion of the partition function 

Solving the transfer equation by means of an inverse temperature power series expansion 



CK{x,y;p) = J2''n{x,y-K)P" 



n=0 

the PDE equation becomes tantamount to the recurrence relation between the densities Cn{x,y; K): 

[nI+ixk-yk){dkI-^Qk)]cn{x,y;K) = [AI-QA-V]cn-i{x,y;K) , n>l 

to be started from: cq{x, y; K) = I. 

The coefficients of the asymptotic expansion for the partition function are obtained through inte- 
gration over the whole plane of the Seeley densities: 



Tre-^^ = ^ E /5"cn(if) , Cn{K) = <fx [Cn]aa{x, x; K) 

n=0 a=l 



Because 

ci{x,x]K) = -V{x) , 
and since indP is independent of /3, we find at the /? = -infinite temperature- limit : 

The recurrence relation gives us the second vortex Seeley density: 

C2(x, X- K) = -\a V{x) + ^Qk{x)Qk{x)V{x) - \dkQk{x)V{x) + \Qk{x)dkV{x) + ]-V''{x) 

The determination of higher-order densities becomes more and more involved. To make the problem 
more tractable we introduce the following notation: 

("--)cf(f) = linr^^^^^^i^P^ , [c„].,(f,f;if) = (°'°)cf(f) . 
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Thus, at the y — > x Umit the recurrence relations between densities and partial derivatives of densities 
can be written in the compact form: 



ai W 02 
r )\t 



^— ^("^-'•+^'"^-*)cf (x) + 



N ai a2 

d=l r=0 t=0 

pir+t r)ad 
_l_ 'J V2 (ax-r,a2-t+l) ndb 

dx\dx\ 



+ 



d=l r=0 t=0 



-•'"»-'>C£i(i) + 



d=l r=0 f=0 

N a2 a 



d=l r=0 t=0 

to be solved starting from 

co{x, x; K) = I 



t ) dx\dx2 
t y dx\dx' 

a 

t J \ r J dx\dx2 



(a,/3)(^afe(^) = 0, if Q ^ 0,and/or/3 / 
(o,o)Qa(-) = 1, a = 1,2, •••,Af 



Knowledge of C2(x, x; requires knowledge of all the densities and their derivatives shown below: 

(2.o)C7o (i'i)C7o (0'2)Co 

(3'0)C7o (2,1)^^ (1,1)C7o 



(3,1) 



Co 



(2,2) 



Co 



(l'3)Co 



(2.0)Ci ('''^Ci (0,2)(^^ 



(0,0)^2 

In general, the number of derivatives and densities required to compute the n*'^-order density is: 

= 16 ■ 1^ ^-^^^-^ ~ = -(n + l)(n + 2)(4n + 3) 

Evaluation of (°'0)Cf''(f) requires knowledge of 4032 local coefficients !!!. The Seeley coefficients are 
then obtained by numerical integration of the Seeley densities over the plane 



r ^ 

Cn{K) = (fxY^ (0,0) ^aa(^ 
a=l 



Note that the upper delta-shaped wing array is fixed by the initial conditions set to start the recurrence 
relations. Thus, strictly one could skip computing the j = n + 1 coefficients in the sum. 
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8.4 The Mellin transform of the asymptotic expansion 

A good approximation to the generalized zeta functions of both K and is given by the Mehin 
transform 



0<{s) 



1 1 1'"^ r 

Jo r(s) Jq 



apphed to the high-temperature expansion of the partition functions 

The generahzed zeta functions are thus divided as sums of meromorphic -high-temperature regime- 
and entire -low temperature regime- functions of s: 

Ck{s) = FfrEE/ d/3r+"-'c„(i^)e-^ + -f- / Tr^e-'^^d/? 



1 _ pi 1 poo 

^ ^ n=Oa=l •^'J V ; --i 

+ ri-l,l] 1 
47rr(s) r(s)' 



, ^.-.lls + n — 1,1] 1 „ , 



n=0 



1 /"l 1 /"OO 



n=0 



We shall neglect the entire parts B{K) and B{K^) and keep a finite number of terms, iVg, in future 
use of these generalized zeta functions for the regularization of ultraviolet divergences. 

8.5 Zeta function regularization 
8.5.1 Tadpole /self- energy graphs 

The contribution to the one-loop vortex mass of the counter-terms induced by the scalar and vector 
fields 

/^M^t =2hml{l) j cPx[{l-\s{x)\'^)-l + l] , AM^t = hml{l) j cP x [^k^k - Vk{x)Vk{x)] 

is proportional to the integral /(I) and, hence, ultraviolet-divergent. To regularize this integral, we 
apply the residue theorem to integration in the complex kQ-jAaxie of /(I) and note that on a square of 
area m?Lp' the integral becomes an infinite sum over discrete momenta: 

1 /■ d^/fc 1 _ 1 1 V- 1 

Thus, this integral is formally the generalized zeta function of the Klein-Gordon operator evaluated 
at s = ^, and we shall take 

^ ' 2m2L2''^o°4^ Svr Til) 47r 
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as the regularized value of /(I). In this way, we find the following contribution to the one- loop vortex 
mass shift: 



Vk{x)) = J d'x [2(1 - \s{x)\') - Vk{x)Vk{x)] 

Contrary to the kink case, which is a one-dimensional problem, a finite answer is obtained in the 
regularized integral via the associated zeta function. The reason is that in this two-dimensional 
problem the physical limit s = ^ is not a pole and only finite renormalizations will be necessary. 
Nevertheless, to keep the procedure as unified as possible we also define the mass renormalization 
corrections as a meromorphic function in the complex s-plane: 

^Mv{s) = ^ (41 Cko{sMs{x), Vk{x)) , AM^ = hm AM^{s) . 



8.5.2 Vortex Casimir energy 

The divergent vortex and vacuum energies 

AEv = —Tt*K2- —Tt* {K'^)2 , A^o = ^Tr " ^Tr {K^) 2 
can be regularized in a similar vein 

^Evis) = ^ [^Jms) - CKois)} , A^o = ^ {Ck.{s) - Cko{s)} 

Recall that 

. , , m^L^ l[s-l,l] , , . m'L^ l[s- 1,1] 

^^^^'^ = — ' ^^0-^^) = ^;^ TisT 

neglecting the entire functions B{K) and B{K^). Thus, 

/i^/^i^yj 21 ,r,os~i , \^ [ rT^\ rT^G\i -f[s + n- 1,1] 



■^0 n=l 



2 \my r(s) Vo ^ ^ ^ 47rr(s) 

where the 21 zero modes have been subtracted: i.e.. 



AMi°)^ = -^l = -0.56419mm 



is the contribution of the 21 zero modes to the one- loop vortex mass shift. 

The zero-point vacuum renormalization, however, amounts to throwing away the contribution of 
the cq{K) and co{K^) coefficients. The physical limit is 



AM{; = lim AM{f{s) , AM^ = lim AAf^(s) 
giving the vortex Casimir energy. 
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8.6 The high-temperature one-loop vortex mass shift formula 

The contribution of the ci coefficients to the vortex Casimir energy is: 



The Seeley densities are respectively scalar, ci{x,x; K*^) 
ci{x, x; K) 



l-|s(^)|' 



2ViS2 
2V2S2 



V 



2ViS2 2V2.S2 %{l-\s{x)\'')~VkVu 
-2V1.S1 



47rr(s) 

s(x)p, and 4 x 4-matrices: 

-2ViSi \ 
-2V2S1 


l{l-\s{xr)-VuVu) I 



-2V2.S1 

Thus, the first Seeley coefficients due to normal and ghost fluctuations are respectively 

tr j (fxci{x,x;K) = j d^x [5(1 - |s(f)p) - 2Vfc(f)Vfc(x)] 

such that 



ci(/^«) 



Am'^^^{-1/2) 



hm 7[-lAl] 
exactly kills the contribution of the mass renormalization counter-terms 



hm 



AMi^{l/2) = — -J:{s,Vk) 
ovr 



7[-iAi] 
r(i/2) 



as expected. In the planar Abelian Higgs model all the particles are massive and we have set our finite 
renormalization prescriptions in such a way that the quantum corrections vanish at the limit where 
the masses go to infinity. 

We finally obtain the high-temperature one-loop vortex mass shift formula: 



AMv 



hm 



-.No S 9/ 

^ n=2 ^ 



The final form is a polynomial expressions in incomplete Gamma functions times the heat-kernel 
expansion coefficients, starting from the second-order coefficients. By cutting the expansion at a finite 
number A'^o we admit an error - besides the rejected entire parts - proportional to j[No — ^, 1] 

for Nq large. 



8.7 Mathematica calculations 

8.7.1 The mass shift of superposed vortices for vorticities 1 = 1, 1 = 2, 1 = 3, l = A 

Plugging the spherically symmetric solutions into the Seeley densities, the coefficients can be calculated 
in a Mathematica environment through numerical integration. In the next Table the result is shown 
for low vorticities: / = 1, 2, 3, 4. 





/ 


= 1 


/ = 


2 


n 


Cn{K) 


c„(i^^) 


Cn{K) 




2 


30.36316 


2.60773 


61.06679 


6.81760 


3 


12.94926 


0.31851 


25.61572 


1.34209 


4 


4.22814 


0.022887 


8.21053 


0.20481 


5 


1.05116 


0.0011928 


2.02107 


0.023714 


6 


0.20094 


0.00008803 


0.40233 


0.002212 
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I 


= 3 


I = 


4 


n 


Cn{K) 




Cn{K) 




2 


90.20440 


11.51035 


118.67540 


16.46895 


3 


36.68235 


2.60898 


46.01141 


4.00762 


4 


11.69979 


0.46721 


14.64761 


0.77193 


5 


2.86756 


0.067279 


3.58906 


0.11747 


6 


0.566227 


0.0079269 


0.667202 


0.01620 



Note that for A^o = 6 the differences between nearest order coefficients are aheady very small, rein- 
forcing the good convergence properties of the high-t expansion. The general formula thus gives the 
one- loop vortex mass shifts, providing the numbers shown in the next Table as a function of A^o- 





AMy(iVo) 

I = 1 


AMy(iVo) 

I = 2 


AMy(iVo) 

I = 3 


AMy(Aro) 
/ = 4 


2 


-1.02951 


-2.03787 


-3.01187 


-3.97025 


3 


-1.08323 


-2.14111 


-3.15680 


-4.14891 


4 


-1.09270 


-2.15913 


-3.18208 


-4.18014 


5 


-1.09427 


-2.16212 


-3.18628 


-4.18534 


6 


-1.09449 


-2.16257 


-3.18690 


-4.18606 



It is also remarkable to realize, as shown in the next Figure, that the mass shift is almost linear in 
l\ i.e. the mass shift for / vortices is almost equal to / times the mass shift for one vortex. 



I 


AMy /hm 


AM 


1 2 3 


1 


-1.09449 


-1 




2 


-2.16257 


-2 




3 


-3.18690 


-3 




4 


-4.18606 


-4 





8.7.2 The mass shift for solutions with two separate vortices 

The same procedure also provides the quantum corrections for two- vortex solutions with intermediate 
separations d = 1, d = 2, and d = 2> between centers, shown in the next Table 





d=l 


d = 2 


d = 3 


n 


Cn{K) Cn(K^) 


Cn{K) Cn{K'^) 


Cn{K) C„(K^) 


2 
3 


61.0518 6.81277 
25.6137 1.33822 


58.3359 6.46609 
24.5050 1.23466 


57.3420 6.03872 
24.1187 1.02031 



No 


AMv{No)/hm 


AMv(No)/hm 


AMv{No)/hm 




d=l 


d = 2 


d = 3 


2 


-2.03770 


-1.99798 


-1.98848 


3 


-2.14095 


-2.09695 


-2.08672 



In this case, the results are less precise for two reasons: First, the non-symmetric approximate solutions 
are less reliable than the superimposed vortex solutions. Second, the first-order equations are not 
enough to solve for field derivatives in terms of the vortex field profiles. Accordingly, we have given 
corrections only up to A'^q = 3. 
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9 APPENDIX I. Zeta functions and Casimir effects 



9.1 The free scalar field in (l+l)-dimensional space-time: quantum vacuum en- 
ergy 

The action for one scalar field (j){y^,y'^) £ Maps(]R-'^'^, M) is: 



either in dimensional = {y^,y^) = ov non-dimensional = {x,t) coordinates. The general 

solution on the interval [—'^, of non-dimensional length mL of the Klein-Gordon equation 

^{x,t)-^{x,t) + (l){x,t) = , (13) 
with periodic boundary conditions: (f>{—'^,t) = is the plane-wave expansion: 

4>{x, t) = J ^ ^ (aneyip[i^/X^t - i^x] + o* exp[-i + i-^x] 

Here, Xn = ^ + I, R = are the eigenvalues of the differential operator Kq = —-^ -|- 1 in the 
space of functions / : ^ C from a circle of radius R to the complex line. The eigenfunctions 

2 

X T~L 

Koexp{in—} = Xnexp{in—} , = — + 1 , neZ 
form a complete orthonormal set: 

ml 

fn{x)0 = —=: • ex.p[i—x] , / dx f*{x)fi{x) = 6nl 
ymL -K -'-^ 

The classical energy is: 

because of the orthonormality relations and the fact that the kinetic and potential energies of terms 
of the form a-nCin or a*_^a*^ cancel (Jeans theorem). 

Canonical quantization leads us to trade Fourier coefficients by operators satisfying the commuta- 
tion rules: 

[aL4]=0 ' [ol.,ai]= ^ni , [an,ai] = 
The quantum Hamiltonian operator becomes: 



The vacuum state is annihilated by all destruction operators, d„|0 >= 0, Vn, and the quantum vacuum 
energy is: 

=< 0|//|0 >= — • 



2 

ne2 



This divergent quantity can be regularized by using the zeta function regularization method and 
expressed in terms of the Epstein zeta function: 

^v(^) = — <kM = — -Y. T^^r—r^ = — ■ E{s, 1|^) . 
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9.2 Inserting plates: quantum Casimir energy 

The Casimir effect measures the quantum vacuum energy when two plates are inserted at x = and 
X = ma,L >> a, with respect to the quantum vacuum energy when the two plates are absent. Thus, 
we must deal with the spectrum of Kq with the Dirichlet boundary conditions: 

4>{ = 0(— ,i) = , (/)(0,t) = (/)(mo,t) = 

The eigenfunctions of Kq for the Dirichlet boundary conditions of the Casimir set-up are of three 
types: 

1. 







mL 
2 



dxf<{x)f<{x)=6ni , neN 



2. 



2 7T7X Tt'^ 

/J5^(x) = \ • sin X ■ "dima — x)'d(x) , A!5^ = — t^—^t • + 1 

ma ma m'^'a'^ 

r dxf<>{x)f<>{x) = 5ni , neN . 
Jo 



3. 



r->/ \ / 2 27rn ^ Att^ 2 

f„{x) = 4 / — • sin — -x ■ r/ix — ma) , X„ = — ttft ttt ■ n +1 

•'"^^ Vm{L-2a) m{L - 2a) ^ ) ^ n rn^{^L-2a)^ 



mL 

dxf>{x)f>{x) = 5ni , neN 







Plane waves in the Casimir set up move in three disconnected regions: 
1. 

' 2h 



')<(X 



t) = \j^l-Yl ("nexp[iyA< t] + {a<yexp[-i^/X< t]^ ■ f<{x) 



1.0 




{x,t) = ^I—.Y^ (a<>exp[iy^ t] + {a<>rexp[-iy^ t]) ■ f< 
ma ^-^ V / 

nGN 



3. 



(/)>(x,t) 



2h 



m(L — 2a) ^-^ 



'a>exp[iyA> t\ + (a>)*exp[-iyA> t\ \ ■ />(x) 
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Accordingly, the Casimir classical energy is: 



Canonical quantization proceeds by requiring the commutation rules: 

[(a<)t,a<] = 6ni , [(a<>)t,a<>] = , [(a>)t,a>] = 

and any other commutator between the creation and annihilation operators equal to zero. The quan- 
tum Hamiltonian for the Casimir set up is therefore: 

The vacuum state is annihilated by all the annihilation operators: 

a<|0(a) >= a<>|0(a) >= a>|0(a) >= , Vn 
and the quantum Casimir energy Ec{a) = Ev{a) — Ey =< 0{a)\H{a)\0{a) > — < 0|ff)|0 > is: 

Ec{a) = ^-(Y^ (^/A< + ^/A<^+^/A>)-E ^) 

VnGN neZ / 



9.3 Zeta function regularization 

We regularize the divergent quantity Ev{a) =< 0{a)\H {a)\0{a) > by means of generalized zeta func- 
tions l|: 



fvm / 

Ev{a, s) = — [CK</Dis) + Ck<>/z)(s) + Ck>/d{s) 



1 



Defining R 



Ck<>/d{s) = 

2 _ m?L'^ ^2 



Ci^>/D(«) = E 



n=i i-m-!,!^ J 

= 'irn^o? , and T?.^ = ^4,^^"^ ? t'^^ three generalized zeta function are given, 



via Mellin transform, by Epstein zeta functions: 



Cxo(s) 

Ci^</D(s) 
Ci^<>/D(s) 

Cx>/d(s) 



1 /"OO °° 2 1 



1 / 1 
2 



1 / 1 



2 r(s) 



1 / 1 
2 



- E\s,l\^] - 1 

2 V ^ ' '7^2J 



By the notation C/c/dCs) we wish to stress the fact that the spectrum of the differential operator K is considered 
with Dirichlet boundary conditions. 
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Application of the Poisson summation formula provides the expressions: 

OO r\ on 





r(s) Jo 


Ck</d{s) 


1 

= 2 




Ck<>/d{s) 




( z 
\2T{s) 


Ck>/d{s) 




[ r(s) 



l=~oo 

OO 



poo ^ 

/ df3p'-2 -e"^ V e 
JO ,_ 



( = — OO 

OO 



/■OO ^ ^2j2 \ 

/ d/?/3'-2 -e-^ ■ V e"~ - 1 
■^0 Z=-oo / 



/ = — OO 



Therefore, 

CK</Dis) + CK>/Dis) - CKois) = ^ 



E_ TV TZ I r — ^ 



e f 



/ = — OO i = — OO 



and the limit when the length of the line goes to infinity is: 

l\m^[CK</Dis) + CK>/Dis)-CKo{s)) = ^■l\m^[n-R]-J^^dpp'-le-^-l 

ma r(s — 



20F r(s 

Addition of 

r(s-i) ^ 2 ^fziY+-2 1 

where 



1 



2 



is the integral representation of the Kelvin functions, provides the result: 

3hm hm ^2 ( 1 

l=\ 

The physical limit s = — i of the renormalized and regularized quantity Ec{a, s) is finite: 



Ecia, s) = Ev{a, s) - Eyis) = + " E (^J " J " ^|-(^^) 



1=1 



Therefore, the one-dimensional Casimir force is: 



„in, X Ec , X hm \ - 1 dK^ , 
The asymptotic behavior of the Kelvin/Bessel function for a large. 



Ki{2mal) ^ /^—•e-^'^™' 
a — > 00 V 4:mal 
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tells us that: 



3hm I m ^ e 3/im , / m 

ErAo-) „ : h\ • > ^ — = : n\ ■ Lis e 

' 4 V 167ra ^ 4 V IGyra 



2ami 

2 ■ 



where Li^[2;] = X]^;^ ^j, |z| < 1, is the Polylogarithm function. 
9.4 Three-dimensional Casimir forces 

The real Casimir effect is measured in a three-dimensional space. As in the one-dimensional toy model, 
the quantum Casimir energy is the difference between the energies of the vacuum state of a free field 
when two plates divide the space in three regions or there are no plates. With no plates and periodic 
boundary conditions chosen to solve the Klein-Gordon equations for one scalar field, the space is a 
three-torus = x §^ x and the quantum vacuum energy is: 

, 711,712,713 GZ , 71 = 71161+71262 + 71363 , 

^ ^ mLi mLi niL^ 

dk ■ dj = Okj , k.j = 1, 2, 3 , Ri = — — , R2 = — — , i?3 = -7, — 

2tt 2tt 2tt 

If two impenetrable plates are located at the two-dimensional sub-spaces {xi,X2,0) and {xi,X2, ma), 
the cylindrical space = x §^ x I is divided into three zones by the plates. We choose pe- 
riodic boundary conditions in the xi,X2 coordinates, and Dirichlet boundary condition as in the 
one-dimensional Casimir effect in the 0:3 coordinate to solve the free field equations. Therefore, the 
quantum vacuum energy in the Casimir device is: 





p2 ' p2 ' p2 ' - ' \ T?'i n^i '2 A • • \ DZ ' p2 ' -p ' 

11^ 1X2 -ftg y it^ 1X2 Z y ri^ JX2 /v- 



3 



rii, 712 G Z , 713 G N, Tlz = ^^^2tt ' whereas the Casimir energy reads: E(P{a) = Ey^{a) — Ey^ . 
The divergences can be regularized in terms of Epstein zeta functions. 



I -n-l -K2 -K3 



hm 

T 



TP\ II ^1 , ^2 63 ei 62 2^31 

It^ 1X2 JX3 Jti JX2 z 



better written, via the Mellin transform, as: 



/■oo 2 2 

/ d/?/5-iexp[-/3(-^ + ^ + l)] 
Jo J^i ^2 
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-» -» ~* 1 /"oo 2 2 2 



Before taking the limit of large L3, it is convenient use the Poisson summation formula to obtain: 
Zi, ^2 G ^) but /i = ^2 = is excluded. Also, since ^3 G Z and /i = Z2 = ^3 = does not count, we have: 



T?\ 1 I ^1 I ^2 , 63 



«GZ3-{0} 





% = 








/•oo 

Jo 




rGZ3-{o} 


6*2 _^ 6*3 J _ 












+ E 


r°° 5 

/ e-^exp[- 
Jo 




rez3-{o} ' 



^^^^ V ^ rez3-{o} ^ 
Thus, the Epstein zeta functions are given in terms of Kelvin/Bessel functions: 

i?[«,i|§ + §] = ^^- (r(.-i)+ E [^'{Rli'i + Rlil)]'^-K^-s{2n^Riq + Riq) 

^ ^ \ zez2-{0} 

i?r 11 ^1 , ^2 , 63 1 
-"-1 ^2 ^3 

TT^RlR2R3 I ^, 3, 



r(.-^)+ [vr2(i?2/2 + ^2^2^^3;32)].(-|).^3_,(2vr^ii2^2 + ^2^2^ 

^ rez3-{o} 



^[->i|^ + ;^ + 4-^3]= (14) 



^ ZGZ3-{0} 
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K2 /C3 

7^ti^li^27^3 / ^, 3 



\ rez3-{o} 
A magic cancelation occurs in the L3 ^ cxd limit: 

r 1 / 11 ^1 , ^2 , 63 ,,6*1 6*2 6*3 



Lr-00 2 \ ^ ' ' i?2 Jl2 ^2J -^''^-^r2 ^2 ^2- 

7rfi2ii22 r(s-|) ^l/2i/22 r(s-|) ma 

- — 7; — - -Ksj ^ ~~Fr~\ — — 

2 L [s) 2 L (s) TT 

exactly cancels with half the first summand in p4p . There are no divergences left in the physical limit 
s = — I and we finally obtain the renormalized, regularized and physical, 3D Casimir energies: 

^qri/ N 3hm ^ ^ r(s — 1) 
S3P(a, s) = —■ 7rR,R2 ■ 

3hm 7ri?ii?2 



r(s) . 

«GZ2-{0} 



rGza-{o} 

+ ^^.^R,R,. Y [^\Rlll + Rill)]--- ■ Ks{27rJRiq + Rjq ) 



«GZ2-{0} 



Isolation of the part depending on a provides the exact result in a very long, L3 ^ 00, cylinder: 
hm , , v-^ 1 



rez3-{o} 

This expression, involving the Kelvin or modified Bessel functions of order two, is not very useful. 
9.5 The infinite-plate area limit 

Taking also the Li — > 00 and L2 — > 00 limits it is possible to obtain a numerical answer. The key idea 
is to perform sums only over and postpone sums in ni and n2 until the end. It is then possible 
to trade the sum of series by the computation of integrals at the limit where the plate areas go to 
infinity. Thus, we have: 

-> — 00 00 22 
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^ ■J ^ ' ni=— OO n2=— OO " 

OO OO 



OO „2p2,2 



Z3 = -00 



^ ' I ni=— oon2=— OO \ r 2 [^—i 



OO OO 



E E l^j^ + ^-E(fr-^V.(*^) 

ni=— 00 722=— OO \ '3=1 



ni=— 00 n2=— 00 



At the L3 ^ 00 limit, the cancelation between infinities explained above takes place and we are left 
with the result: 



Ef[a,s) = — -{ 



00 00 



E E -Et^r^^^i-.i-fap) 



ni=— 00 712=— 00 yP^^^^^P ^ ^3=1 

which at the physical limit s = —\ is: 



2 



hm 



00 00 



7ii=— 00 n2=— 00 



2 °° 1 S 

-E r •^i(2«"i^3^') + ^ 



vr Z3 

(3 = 1 ^ 



At the limit of very large areas of the plates the discrete momenta become continuous and the series 
become integrals according to very well known prescriptions: 



00 00 



m?LiL2 — > 00 47r2 



00 /"OO 



7ii=— 00 n2=— 00 

Therefore, the part of the Casimir energy depending on a is: 



00 J —00 



dpi / dp2 



hm m?'LiL2 2 



2 47r2 



- • / dpi dp2P-^ -■ Ki{2amhP) 

J — OO J —OO h = l 



Because dpi (ip2 = 27r pdp = 27r PdP, p^ = Pi + P2, the integral of the modified 
Bessel function can be performed to find. 



EW(a) 



h L1L2 1 h 



^ '"^^ 87r2 a3 ^ /4 ^3 

(3=1 3 

and the Casimir three-dimensional force is: 



h L1L2 ^^^^ _ -^^1^2 



720 a'^ 



da 



240 a* 



78 



10 APPENDIX II. Kinks: d=l, N=l 



Kq = —-^ + 4 acts on functions / : ^ C from a circle of radius R = (PBC), to the complex 

plane. The spectral resolution of Kq is: 



Koexp{m^} = A„,exp{m-^} , = ^ + 4 , nS 



and the vacuum energy on the circle is: 



.21 4 



Regularization of this divergent quantity by means of the generalized zeta function leads to: 



2 ' V"i2 



s G 



Cxois) = Tr 



+ 4 



(n2 + 4i?2)« 4 



1. (n2 + 4i?2) 

n=l ^ ' 



10.1 The kink generalized zeta function versus Riemann zeta functions 

To express the generalized zeta function in terms of ordinary Riemann zeta functions one can use the 
binomial series: 



n=l ^ ' ' 1=0 ^ ^ n=\ 



10.2 The kink generalized zeta function versus Epstein zeta functions 

Alternatively, the generalized zeta function is an Epstein zeta function, 



oo 



which via the Mellin transform 



1 °^ fCO 

E{sA\A)=—-- / 

\ ' r, — — < , •'0 



-/3(An2+4) 



and use of the Poisson summation formula 



oo j oo 2;2 

E -^*-^ = ,/^-E »p{-^) 



i=— oo 



A(3 



79 



reads: 

1 r— f /-oo /-oo 2/2 ) 

E{s,4\A) = ^y^j-lJ^ dpp'-2e-''' + 2j2j^ d/S^-^ exp{-4/3 - -^}| 
When i? oo, A = ^ 0, the result of Section §3 is obtained: 

In sum, 

1. The result is the same when we work with periodic boundary conditions on a finite interval and 
allow the length to go to infinity as when physicist's techniques are used to cope with continuous 
spectra. 

2. In both approaches there are contributions proportional to the volume: infrared divergences 
which must somehow be renormalized. 

10.3 The high-temperature expansion of the kink heat equation kernel 

Let us now consider the differential operator 

K = -— + A-V{x) , V{x) = V{x + 2TrR) 
acting on functions from a circle of radius R to C. If V{x) = 0, we find the -fCo-heat kernel: 

oo 2 

KKo{x,y;fi) = ^ exp{-/3(-^ + 4)} • exp{i— (x - y)} 

n=— oo 

= e-^/^exp{-^:^} . f; exp{-A[„ + ,|(y_,)]2} . 

n=— oo 

Using again the Poisson summation formula, 



^ I — 2/2 a T> 

5: e-(-^^)^ = J^. 5: exp{-^-2vrzM ; t = A , v = ^§-{y - x) 

n=—co ' i=— oo 

the i^o-heat kernel can be written as: 



4fl ttR"^ , (y — x)"^ ^ , ttRIIttRI — (y — x)] . 
KK,{x,y,l3)=e-^^ ■ ■ exp{-LA} . ^ exp{ ^ ^} 

For f3 < 1, (high-temperature), we obtain the asymptotic formula. 



and from this, the asymptotic high-temperature expansion for the ii'-heat kernel is derived: 
KK{x,y;P)c^e~^P ' VT" ' """P^"^^!^^ • c„(x, y) /3" . 



n=0 
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10.4 Kink Seeley densities 

We now describe the iterative procedure that gives the coefficients Cn{x,x;K) used in the text. For 
an interesting interpretation of these coefficients as invariants of the Korteweg-de Vries equation, see 
the book on Quantum Mechanics by Perelomov and Zeldovich. 

The recurrence relation 

(n + 1) c„+i (x, y) + ix- y)^^^^^^ + V{xMx, y) = ^^^^ (16) 

comes from phigging the high-temperature expansion into the transfer equation. In order to take the 
limit y — > X properly, we introduce the notation 

= hm ^^^^^ , 

and, after differentiating (jl6p k times, we find 

n + k I ^ \j ) ox^ 

From this equation and ^^'^Cq{x) = limj^^^: = 5^^ , all the ^^^^Cnix) can be generated recursively . 
We finally obtain a well-defined recurrence relation 

Cn+i(x,x) = ^ ^"^^Cnix) -V{x)Cn{x,x) 

n -|- 1 L 

suitable for our purposes. 

We give the explicit expressions of the first eight c„(x,a;) kink coefficients. The abbreviated 
notation is Uk = ^{x), = (^(x) 

ci(x,x) = Uq 

C2{X,X) = -UQ + -U2 

2 b 

, ^ 1 3 1 1 2 1 

C3(X, X) = -Uq + -U2U0 + — Ui + —U4 

, N 1 4 1 2 1 2 1 1 2 1 1 

C4(X, X) = —Uq + ■:^U2Uq + r^U^Uo + -^U^Uq + —U2 + -^UiU-i + —Uq 

, N 1 5 1 3 1 2 2 1 2 1 2 1 1 11 2 

5^ , ) 120 ° 36 ° 24 ^ ° 120 ° 40 2 30 840 360 ^ 

23 2 19 1 1 

+ H U2Ua H U\U^ H lis 

5040 ^ 2520 280 15120 

cel^;, x) = —Mo + Y44"2^io + T^^i^o + ^^^^^o + ^"a^^o + ^^i^s^o + —u^U2Uq + ^^^i^^s^^o 

1 4 1 61 , 43 23 2 5 2 19 

+ 288^^ + 15120"«"° + 15120^2 + 2520^^"^"^ + SOio"""^ + lOOs"^"^ + 252o"°"^^^ 

23 2 19 1 2 11 1 1 

+ 30M0"^ + Tsm^^^^ + 1680^°"^ + 15m"^"^ + 3780^^"^ + 33^^^° 

C7(X, X) = ^^Ul + ^^U2UI + ^^uW, + ^^Ulul + -^^U,U,ul + ^^^2^^ + ^^^5^^ 
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+ 
+ 
+ 
+ 
+ 



1 4 61 o 43 So 

288 ^ 15120 2 2520 ^ ^ ^ 'J ^ i008 ^ 

19 2 2 2 ^ S 

——UoUAUn + — — -UfiUn + —TTTrrrU.Un + -—UtU?, 



— - — ujuAUn H UinUn H ^ — U^Un 

1008 1 'J ^ 332640 10080 ^ ° 

1 



zoo iOiZU ZOZU ±UUO OOZU'iU lUUOU 

19 2 1 . 83 2 2 1 s 31 2 1 

5040 ^ * ^ 5040 ° 10080 ^ ^ 252 ^ ^ 10080 ^ ^ 280 

5 2 23 2 1 19 71 2 1 2 

2016 ^ 30240 u 4 -r ^ ^ ^ ^ 15120 'J ^ ^ 665280 ^ 2016 ^ 

11 61 1 19 1 2 17 

15120 " ^ ^ ^ 332640 3780 " ^ ^ ^ 1663! 

1 1 

U\Ua H U\o 

66528 8648640 



1440^°"^ 



x) 



+ 
+ 
+ 



40320 ° 960 2 ° 720 ° 576 ^ ° 252 ^ 280 ^ 420 ^ ^ ^ " 

31 2 5 2 1 2 1 23 2 2 19 
uoUnUnA unUAUc + ^^^i^e^o + gg4gg4Q ^i2no + ■ 

,2 11 



11 



23 2 2 — 
60480 ^ ° 30240 



'0 



1 fi 17 4 83 2 2 61 2 s 1261 4 43 n 

UnUo H U-1U2 H UnUlUn -\ UnUn -\ Un -\ UnUTUoU' 

4320 ° 8640 ^ 10080 ^ 2 ^ 30240 ° ^ 1814400 ^ 5040 " ^ ^ ' 

227 2 23 o 2 659 2 2 5 2 2 19 s 527 

37800 ^ 2 ^ ^ 3Q240 3 "t 392400 ^ ^ ^ 



+ 
+ 
+ 
+ 
+ 



1261 4 43 2 

71 ^2 "I ^^^o''^1^2^3 

'iozu oo^u luuou ouz'iu 1814400 5040 

227 2 23 3 2 659 2 2 2 2 19 3 527 

37800 "^"^"^ + 30240^0^3 + 302400"^"^ 2016^°"^"" ^ 15120 ^ 15120C 
7939 2 6353 2 1 3 17 s 13 306' 

9979200 ^ 9979200 ^ 1680 ° 30240 ^ 12320 49892 

71 2 1 4 3001 2 13 61 

665280"°"^ + 20160^0^^ + 99^^2-6 + ^^1-3-6 + ^77-0^4^6 



151200 ^ 
13 3067 

UoU-^Uk. -\ U1UAU5 

12320 4989300 

71 2 1 4 3001 2 13 61 

665280"°"^ + 20160^°^^ + 9979200"^"^ + 29700^^"=^"^ + 33264o"°"^^^ 

3433 2 109 19 1501 71 2 

259459200 ^ 498960 166320 64864800 1995840 ^ 

17 2003 1 5 1 2 

33^"°"^^« + 129729600 + 66^8^°^^^^ + 648648 "^"^ + 6™"°^^° 

73 1 1 

U2U^() H u^um H U-IA 

25945920 1441440 259459200 



11 APPENDIX III. Two-component kinks: d=l, N=2 

The operator 



+ 4 







acts on functions / : ^ C © C from a circle of radius R = (PBC), to complex isospinors. The 
spectral resolution of Kq is: 



exp{4 (exp{^4^.x} 



^(2) 



(n(2)) 2 



+ 4 



n 



(2) GZ 



and the vacuum energy on the circle is: 
hm 1 







2 R 



i(i)=-oo 



(l))2+4^2 



00 

E [( 



+ > \[n 



(2))2 + ^2^2 
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Regularization of this divergent quantity by means of the generahzed zeta function affords: 



AEois) 



2 \m 



2 ; • ^ • Cko{s) 



see 



Tr 



+ 4 



E 

nt-'-) — oo 



((n(i))2 + 4i?2)^ 



2R 



2s 



E 

i(i)=i 



l(2)= — oo 

1 



((n(2))2 + a2i22). 



((n(l))2+4/?2)s Cj2s 



l(2)=l 



((n(2))2 + a2i?2). 



11.1 The TK2 kink generalized zeta function and Epstein zeta functions 

The generalized zeta function is the sum of Epstein zeta functions, 

which via the Melhn transform, 

1 POO 

E{s,4\A) = —-. ^ / 



'/3(An2+4) 



and use of the Poisson summation formula 



E 



i=— oo 



A/3 



reads: 



r(s) 



2/2 



When R ^ oo, A = ^ 0, the result of Section §4 is obtained: 



E{s,4\0) + E{s,a'^\0) = lim 



1 



A-*o r(s) V -4 



A/3 



lim 



mL , 1 



- + 



d/3/3"-2 [e-^^ + e"" ^ 

1 . r(s-i) 



CJ 



2s-lJ 



r{s) 



11.2 The high-temperature expansion of the TK2 kink heat equation kernel 

We now consider the differential operator 



K = Ko + Vix) 



dx'^ 



+ 4- T/"(x) 








dx^ 



+ a'-V^^(x 



V{x) = V{x + 2ttR) 
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acting on functions from a circle of radius to C ® C. If V{x) = 0, we find the i^o-heat kernel: 
( Eni^^^-o. exp{~Pi^ + 4)} . cxp{z41(x - y)} 

V Eni^^=-oo cxp{-/3( + a^)} ■ exp{*41(x - y)} 

Using again the Poisson summation formula, 

CO j — °° 2/2 o r> 

n=— oo i=— CO 

the -f^o-heat kernel can be written as: 



KKo{x,y]P) = • exp{ — — } X 

-4/3 y-°° I- 7rfi;<^[7rfl;<^-(jy-x)] - 



^ e Z^;(2)=_oo sxpi g i 



For /3 < 1, (high-temperature), we obtain the asymptotic formula: 



^i^o(x,l/;/3) - ^-^ • exp{-^^} • ( ^ j • (l + 0(exp{--})j 

and from this, the asymptotic high-temperature expansion for the i^T-heat kernel is derived: 



~ • exp{ • |^ ^ ^_^2^ j • 2_, c„(x, y; iT) /3 



11.3 Conserved charges of the N x N matrix KdV equation 

We now describe the iterative procedure that allows us to compute the coefficients c^^{x,x; K{c)) 
used in the text. We shall present the diagonal coefficients as the invariant charges oi a N x N matrix 
generalization of the Korteweg-de Vries equation. In this subsection we shall derive the formulas for 
arbitrary dimension N of the target space. We shall also deal with a slightly more general situation 
in which all the particles have different masses, i.e. the vacuum fluctuation operator is of the form: 

\ 


V ••• -£+vf^^ ) 

The kink moduli space for N scalar fields is of dimension N; by c we shall denote collectively the 
A'^ — 1 integration constants that determine the kink orbits. 



+ V 



(1) 



(2) 
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We start from the recurrence relations 

dx 



{n + l)ci^,{x, y; K{c)) + {x - y)^{x, y; K{c)) 



(x, y- K{c)) + {vf^^ - vfs))c^''{x, y; K{c)) 



dx"^ 

N 

- Y.V^^{x)c^^{x,y;K{c)) 
C=l 



coming from plugging in the high-temperature expansion 



^ 00 

C^^{x, y) = -j= . J2 y, K{c)W 



n=0 



in the transfer equation: 



d X — y d d"^ 2 2 



. N 

\ ■ C^^{x,y) + V^^{x)C^^{x,y) = 0^^ 
^ c=i 



In order to take the limit y ^ x properly, we introduce the notation 



ptk AB 

«C„^^(x) = lim^(.,,;K) 



^^)C^^{x) = hm = 5^H^^ 

y^x Ox'' 



and, after differentiating (I17p k times, we find: 
1 



^'^^^C^^{x) + {vl^-vl^^)^^)C^^{x) 



'+^^~' n + k + 1 
The recurrence relations become 



TV 



when ?/ — > X. From this equation and knowledge of ^^^Cq^{x), all the Seeley coefficients can be 
computed recursively. For instance, the lowest-order Seeley densities are: 



c^^(x,x;i^(c)) = 5 



AB 



ct^ix,x;K{c)) = -y^^(x) 



AB/ 



4''{x,x;K{c)) 



l^yAB 
6 dx"^ 



4^{x,x;K{c)) 



1 dV^^ , . 1 
(x 



c=i 

N 



60 dx^ 



c=i ^ 



TV 



C=l 



C=l D=l 



1 



+ T7;iV(A)-V(B)) 



The diagonal terms c^^{x,x; K) are the densities giving the infinite conserved charges of a x 
matrix Korteweg-de Vries equation, namely: 



dV . 
'dt 



dV 



dV 



ix,t) - 3 ( V{x,t) — ix,t) + —{x,t)Vix,t) ] + -T^ix,t) = 



dx 



dx^ 



(17) 
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where the matrix potential now evolves in "time " t, V = V{x,t). The reason is that this equation 
can be written in the form: 



Lt+[L,M]=0 



^2 



with B{t) an arbitrary matrix of functions of time. Therefore, standard arguments guarantee that the 
time evolution ruled by (I17p produces a uniparametric isospectral transformation of the Schrodinger 
operator L. Because the integrals c^^{K) are determined by the spectrum of L, their invariance 
follows. 

11.4 TK2 kink Seeley densities 

From the transfer equation for C{x, y; (3) in a scalar field theory with two phonon branches of gaps 4 
and 



^ + ^ ^ - ^ + V^" (-) - 4 ) Cj,^ (X, ^) + (X, y ; /3) 



d 




y d 


d^ 


d(3 


p 


dx 


dx"^ 


d 


X — 


y d 






p 


dx 


dx"^ 


d 


X — 


y d 


52 


dp ^ 


p 


dx 


dx"^ 



rll I 



<12t 



^22/ 



rV2(„\rf\2( 



we derive the recurrence relations for the Seeley densities: 

(n + l)ciVi y; K) + (X - y) - ^^^^^ + ( - 4)0^1 (x, y; K) + V'^c^^ {x, y; K) = 



dx dx 



(n + l)c„+i (a;, y; K) + [x - y) 

(n + l)c22+i(a;,y;Jf) + (a; - y)^^^ 
Denoting 



a. + - <^ ^) + ^^"c-(x, J/; ) = 



('=)(:7l^'(x) = lim 



^''^i{x,y■,K) 



(18) 
(19) 
(20) 
(21) 



differentiating the recurrence relations above k times, and taking the y ^ x limit, we obtain new 
recurrence relations between the derivatives of the Seeley densities when y = x: 



ciUx)=^{^''^'^c:\x)-Y: 



'^n.k I 



k 



(k) 



Cn + lix) — 



3 = 

('=+^>C,?(.)-^^ 



1 



(fe) 



Cn+i(2;) 



<'n,k I 



(fc+2)^22 



cr(-)-E 

J=0 













){k 










dx3 






0^ ){k 


dxi 



'C7"(x) + 



'C„ (x) + 



)cr(x)+ 



9x3 



9x3 



9x3 



C„ (x) 



where m„ ^ = n+A;+l. Moreover, the infinite temperature condition requires that: ('=)Co'(x) =6'''^6'^. 
Thus, the recurrence relations become: 

1 



.11 , 

-n+l 



{x,x) 



n + l 



(2)c"(x) - {V'Hx)-A)ci\x,x)-V'\x)cf{x,x) 
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1 



-n+1 



„22 



n + 1 
1 

n + 1 
1 



^'^Cl\x) - {V'\x) - A)cf{x,x) - V'H^)ci\x,x) 



^'^Cl\x) - {V'\x) - a')cf{x,x) - V'\x)ci\x,x) 



which must be solved iteratively. Up to third order, the Seeley densities are: 



c"(x,x) = 
cf{x,x) = 



-V'^x) 
-V'^ix) 
-{V'\x)-a') 



c"(x,a;) 
cl^{x,x) 

^2 (^i^) 



2ta11 



6 dx^ 
1 a2yi2 



1, 



6 dx^ +-,V''ix)[V'\x)+V^\x)~2a^] 



„ii 



„22 



60 dx^ ^6^ ^""^^ ^ +6*'^ dx^ ^12 ax 



1 /^T/ia 1 



1 



(Fi2(.^))2(^ii(^)_4) 



1 9^22 1 



Q2T/22 1 ayl2 



1 91/22 gy22 

12 i9a:; 5x 



-liV''ix)nv'\x)-a') 



12 APPENDIX IV. Self-dual Vortices: d=2, N=4 



At the self-dual limit = 1 there are two PD operators ruling the small fluctuations around the 
vacuum of the Higgs, Goldstone, vector, and ghost fields. In the ii-gauge these operators are: 



a 









+ 1 
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d 



dx\ dx\ 



+ 1 





, 92 _ a 






+ 1 






92 d 



dxi 



dxt. 



+ 1 
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d 



dx\ 8x2 



+ 1 
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Kq acts on isospinors f : x — > from a torus of area A = AtiE? = m'^L'^ to C^, whereas Kq 
acts on scalar functions f : x — > C from the same torus to C. The eigenfunctions of both Kq 
and Kq are plane waves of discrete momenta: 

T^AA I •'''1 •'''2 1/4 \ ) • '^1 ■''2 I A 

° exp < + ^^-^2 }u = A^(A)exp <^ + ^^-^2 ? u , 

Kq exp l^-^a^i + «-^2;2| = A^exp Y^^i + ^"^^2]- 

The eigenvalues are 

and the vacuum energy on the finite torus (no infrared divergences) is: 

Ai?o = — -Tr/^o^- — •Tr(Ko°)2 =3— E + # + ' 

ni=— 00 n2=— 00 

coming from small fluctuations on the vacuum of the Higgs field and the two physical polarizations of 
the vector field. The fluctuations of the temporal polarization of the vector field and the Goldstone 
field are canceled by the ghost fluctuations, thus restoring the unitarity, that was lost in the combined 
Weyl/ii-gauge. 

Zeta function regularization of this ultraviolet divergent quantity leads us to replace ASq by the 
meromorphic function of the complex variable s: 

^^°(^)=K&) ■^■^^'o"^'^='K^) '^'^w^^^s 

and take as the regularized finite value of the vacuum energy the value obtained by analytic continu- 
ation of the zeta functions. 

12.1 The vortex generalized zeta function versus Epstein zeta functions 

Therefore, the vacuum energy is regularized by means of Epstein zeta functions: 

4 ^ 4 ^ ^ 

Ciro(s) = E E -ZiATjA) — -ariA) = E + 



Via the Mellin transform 

and use of the Poisson summation formula 



2^ e ^"15^ = — - • 2^ e /3 , / = Ziei + Z2e2 , /i,/2GZ 
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the Epstein zeta function reads: 



-2,-/3 



+2 E 



d(5(5' 



At the hmit of infinite area only the first term survives and we obtain: 

m^L^ r(s - 1) 



E(g,l|0) = lim R^T^ ^^l, , 

^ ' ' ^ R^oo r(s) L^oo 47r 



hm 



T{s) 



Thus, in the Euchdean plane the regularized vacuum energy is: 



,2 \ * 



lim 3— — 7T ■ U ■ lim 



m 



L^oo 47r 



■L2 r(3- 1) 



Note that 



-| is a regular value of AEq{s). 



12.2 The high-temperature expansion of the vortex heat equation kernel 

Let us consider now the PD differential operators K = KQ + Qk{x)-^ + ^(2;) and K'^ = Kq +V^{ 



Qk{xi,X2) = Qk{xi + 2'kR, X2 + 2ttR) 



/ 



-2Vk{x)^ 

\ -2Vk{x)^ 



x): 
\ 

J 



( w 



V{x) 



1 











\S\X] 



1 



-2ViS2(x) 

-2V2S2(X) 



-2ViS2(x) -2V2S2(f) i(|s(f)|2-l) + yfc(x)T4(f) 



2ViSi(x) 
2V2Si(x) 


\ 2ViSi(x) 2V2Si(x) |(|s(f)|2-l) + T4(f)I4(x) / 

y(xi,x2) = F(a;i+27ri?,X2+27rii) ; = |s(f)p-l , ^'^(xi, X2) = F'^(xi+27ri?, X2+27ri2) 

K acts on isospinors f : x — > from a torus of area A = 47ri?2 = m2_j,2 whereas K"^ 

acts on scalar functions f : x — > C from the same torus to C. The Kq and Kq heat kernels 
are: 



K^t{x,mP) 



e ^exp{- 



4/3 



(x - y)} 



}• E exp{-A.(^(^)+,^(f_y-)).(^(A)+i^(^_y-;))} 



KKi^{x,y,(3) 



e ''expl- 



ii2 

\x-y\'^ 
4/3 



E exp{-/3(^ + 1)} •exp{i^(f - f)} 



} • ^ exp{- 



Using again the Poisson summation formula, 



E 



-t\n+v\'^ = - X] 



exp 



t2 



i? R 

(n + i^ix-y)) ■ {n + i^{x - y))} 
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the Kq and Kq heat kernels can be written as: 



-e ^exp 



/5 



4/3 



^ exp 



e ^exp 



4/3 



vrffi • {ttRI -{y-x)) 



For < 1, (high-temperature), we obtain the asymptotic formulas 



KKoix,y;f3) -- 
KKG{x,y;P) 



-e ^exp 



ttR^ 



4/3 



I + 0(exp(--)I) 



^e-^exp{-^}.(l + C.(exp(-^)) 



and from these, the asymptotic high-temperature expansions for the K and heat kernels are 
derived: 



Ki^^(f,y;/3) = 



-e '^exp 



ttR^ 



e ^exp 



4/3 
4/3 



n=0 
2 ^ °o 



^c;^^(x,y;K)/3" 

1=0 

oo 

Y,c^{^,y;KW 



n=0 



12.3 Spherically symmetric vortex Seeley densities 

Defining 9i = {I — 1)6 and A(r) = [1 — the Hessians K and for spherically symmetric vortex 

solutions read: 



K 



( K^^ ^13 ^14 

^21 ^22 ^23 ^24 

^31 ^32 ^33 ^34 

V A'41 ^42 ^43 ^44 J 

52 1 a2 



K 



G 



dr 



1^ 

r2 



^11 = k22 



dr 



.2 5^2 



^24 ^ ^42 



A(r)sin6l« = -K^^ = -K 



-31 



K^^ = K 



= ^21 ^ Q 



41 ^ 2^A(r)cos^, = -^23 



-K 



32 



^33 ^ ^44 



92 



1 



5r2 r2 502 2 



/2a(r) 1 



^34 ^ ^43 



' r2 50 



where the first-order equations have been used to write field derivatives in terms of field profiles 



ds^ _ Ifjr) 
dxi r 

dsi ^ Ifjr) 
8x2 r 

ds2 _ Ifjr) 
dxi 

dxi 



r 

I cos 26 



cos6 cosl6{l — a{r)) + sin sin /0] 
sin0cos/6'(l — a(r)) — cos sin 
cos 9 sin W{1 — a{r)) — sin 6 cos W] 
lcos'^6{f{r) - 1) 



sm t/ cos ( 



r2 2 



ds2^ _ Ifjr) 
dx2 r 
dVi 
dxi 
dVi 

8X2 

8V2 

8x2 



sin0sin/0(l — a{r)) + cos cos 



-I COS 26^ + ]- sin2 6{f{r) - 1) 



sm a cos ( 



2//(r)a(r) 1,^2 



+ 2(/ 
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The Seeley densities for spherically symmetric vortex and ghosts are: 



tr[ci](x,x;K) = 5[l-/2(r)]-^/V(r) 

tT[c2]{x,x;K) = Y^{37r^ + 4/V(r) + 8(7/^2 -8r^)/2(r) + 27rV^(0- 

- 8lr^a{r)[-l + (1 + 131) f{r)] + 8l'^a'^{r){-2 - 3r^ + 9r^f{r))} 

tr[c3](x, x; K) = {-Afa^{r) - 41^r'^a\r)[U + (-132 + 167/)/2(r)] + 4/^a^(r)(20 + Or^ + 32r2/^(r)) 

- 2/r2a(r)[-4(16 + Or^) + +(64 + 96/ - 472/^ + 344/^ + 88/^ + 243lr^)f{r) 
+ (-52 + lOdiy^f^ir)] + fa'^{r)[-256 - UAr"^ - 117r^ 

+ 2r2(88 - 548/ + 516/^ + 183r2)/2(r) + 99r^/^(r)] + r'^[r'^{-16 + ISlr^) 
+ (-320/2 + 160/^ + 32r2 + 48/r2 - 321r'' + 8/^(20 + 39r'^))f{r) 
+ r2(-16 - 48/ + 44/2 ^ 199^2^^.4^^^ _ 29r^f{r)]} 

ci(x,f;ir^) = l-f{r) 
1 

C3(x,x;K^) = ^-rjlOr''- [-32/3 + 16/'' + 8/r2 + 23r^ + 16/2(1 + r2) 
60r* L 

- 8/(-12/2 + 8/3 + + 4/(1 + r^))a{r) + 16/2(1 - 6/ + 6P + r2)Q2(r) 
+ 32(1 - 2l)l^a\r) + 16/^a^(r)]/2(r) + r2[8/ + 16/^ + 17r2 + 16/2a2(r) - 

-8/(1 + 4/)a(r) ]/^(r) - 4rV^(r)} 

The Seeley coefficients are obtained from numerical integration over the whole plane of the above 
densities when the field profiles /(r) and a{r) of the vortex solutions are plugged in. 



C2ix,x;K^) = --^{[4/2 + 5r2-8/2a(r)+4/2a2(r)]/2(r)-3r2-2r2/4(r)} 



12.4 Renormalization of one-loop divergent graphs in the planar Abelian Higgs 
model 

In this last Appendix we shall present detailed calculations of Feynman amplitudes for one-loop di- 
vergent graphs in the Abelian Higgs model. The results obtained here have been used in sub-section 
§• 7.3. 

12.4.1 The Higgs boson tadpole 

• Application of the Feynman rules gives the following Feynman amplitude to the contribution of 
a Higgs loop to the Higgs tadpole: 

3.2 f d?k i ■ 2 T ( 2\ 

2 J (27r)3 -K? + ie 2 ^ ' 

A combinatorial factor for this graph of ^ has been taken into account. The Mc Laurin expansion 

— = V(-l)"^ 

of 




^ n=0 ^ 



1 + ^ \k^ — ir 



k"^ — K? + ie k"^ — 1 + ie ^ [A:^ — 1 + ie 

n=0 
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= /(I) + f If. /(g, 



allows us to write: 

[K- - ir ■ I ■ 

n=l 

Because all the integrals in the sum - except /(I) - are convergent, we can safely conclude: 

I(k^) = /(I) + finite part 
which is a useful relation if a minimal subtraction scheme is to be used. 

Similar calculations show that the Feynman amplitude of the Higgs tadpole due to a Goldstone 
loop is: 

/""""\ _ 1.2 f ^^fc i 1.2 rn\ 

f — IK, ■ / -. r^TTi = - J 1 

\, 2 J (27r)3 -l + ie 2 ^ ^ 

There is also a combinatorial factor of ^ and the two differences are that the trivalent Higgs- 
Goldstone-Goldstone vertex does not have a factor of three and the Goldstone propagator has 

poles at ko = ±Vkk + 1, contrary to the Higgs propagator with poles at fco = ±\/kk + K? . 
The Feynman amplitude of the Higgs tadpole due to a ghost loop is: 

There is no combinatorial factor and a minus sign is included, as corresponds to all fermionic 
loops. 

The Feynman amplitude of the Higgs tadpole due to a vector boson loop is: 

2'^ J (27r)3A;2-l + ie- ^^^> 

There is a combinatorial factor of ^ and g^'^g^i^i = 3 in (2 + l)-dimensional Minkowski space-time. 

Adding the four summands, the result used in subsection §. 7.3 is obtained: — 2z(k^ + 1) • /(I) + 
finite part. 

12.4.2 The Higgs boson self-energy 

• Identical calculations to those performed in the previous subsection give the Feynman amplitude 
for the self-energy of the Higgs boson due to a Higgs loop: 

= IK'^ ■ / -TT-T^ 7i = IK^ ■ I{k^) 




2 J {2TTf k"^ - K? + ie 



There is also a ^ combinatorial factor and now a four-valent vertex of four Higgs particles replaces 
the former three-valent vertex. 

The Feynman amplitude for the Higgs self-energy produced by Goldstone loops is almost the 
same: 



1 2 f '^'^^ ^ 1 2 / \ 

-IK ■ / r^TT^ = IK ■1(1) 

2 J (27r)3 fc2 _ 1 + ie 2 ^ ' 



The only differences come from a different four-valent vertex of two Higgs and two Goldstone 
particles and a different propagator: Goldstone instead Higgs. 
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A vector boson loop in the Higgs self-energy is easily dealt with. The Feynman amplitude reads: 

ha^"' ■ I '^^"^ - -3? • 

2^ J (27r)3fc2-l + k- ^^^> 

The computation of the Feynman amplitude of the Higgs self-energy due to the only divergent 
two- vertex graph with two external Higgs legs is more difficult: 

\_y ^ J W)^' [{p-k)^ -1 + ieW -1 + ie] 

■ f d^k i{p + k)^ ^. 

V (27r)3[(p-A;)2-l + ie][fc2-l + ie] ' ^^' ^ 

Note, however, that: 

/(p;l)=/(l) + hip) + hip; I) 

d?k Uipk) 



II) + 



+ 



(27r)3 Hp - A;)2 - 1 + ie][k^ - I + ie] 
d?k i(l — ie) 

(27r)3 Hp - A;)2 - 1 + ie][fc2 - I + ie] 



hip J 1) is convergent but it is not obvious that hip) is also convergent 0. We use the Feynman 
parametrization 

1 dx 



ah Jq [ax + 6(1 — x)Y 

to write hip) in the form: 



(27r)3 Jq [A;2 - 2xipk) + xp^ - I + ief 
Changing variables to q = k — xp, 



hip) = dx j 



d^q \ Aiipq) Aixp^ ^ \ -z . 

+ r^r- — a . . in , A* =xil-x)p -1 



(27r 



|3 



[g2 -I- ^2 _|_ ^£j2 Jg,2 _|_ ^2 _|_ ^gj2 

one sees that /2(p) is indeed convergent. Therefore, 

lip; 1) = /(I) + finite part 

Adding the four summands, the result used in subsection §. 7.3 is obtained: — 2(K2+l)/(l)+finite part. 

12.4.3 The Goldstone boson self-energy 

• The Feynman amplitude of the Goldstone boson self-energy caused by one loop of Higgs particles 



is: 




1.2 [ d^k i _ 1.2 J, 2^ 



■^iK^ ■ /(I) + finite part 



'^It could be guessed that l2{p) is indeed convergent because the integrand is odd in the loop momenta k. 
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The Feynman amplitude of Goldstone boson self-energy from one loop of Goldstone particles is: 



in ■ / — 7^—^ = — IK ■ 1(1) 

2 J (27r)3 fc2 _ 1 + ie 2 ^ ' 




The Feynman amplitude of Goldstone boson self-energy from one loop of vector boson particles 
is: 

The Feynman amplitude of the two-vertex Goldstone boson self-energy graph is: 

" J (27r)3 ■ [{p - A:)2 - 1 + ie] [k^ - + ie] 

. f d^k i{p + kf ^ 2 

V (2vr)3[fc2_K2 + ie][(p-A:)2-l + ie] ^MP,^,-^; 

As in the fourth item of the previous subsubsection §. 12.4.2, the divergent integral is the sum 
of three integrals: 

I(p;1)=/(k2) + l2{p;^^)+h{p;K^,l) 

(fk Ai{pk) 



+ 



(27r)3 [{p - A:)2 - 1 + ie][k^ - + ie] 
d?k i{l - ie) 

(27r)3 [{p - A:)2 - 1 + ie] [A;2 - k? + ie] 



Again l2,{p', k^, 1) is convergent and routine Feynman parametrization shows that /2(p; is also 
finite. Therefore, 

/(p; K^, 1) = /(I) + finite part 
Adding the four summands the result used in subsection § . 7.3 is obtained: — 2(K^-|-l)/(l)-|-finite part. 

12.4.4 The vector boson self-energy 

• In the computation of the Feynman amplitude for the vector boson self-energy graph with one 
Higgs loop there is a combinatorial factor of | and we include a factor of three to take into 
account the three possible polarizations of the ingoing and outgoing vector bosons: 

= Sigi^" ■ /(I) + finite part 

The Feynman amplitude for the vector boson self-energy graph with one Goldstone loop is 
identical except that the Goldstone propagator replaces the Higgs propagator: 

^AAAAAAAA/ ^ ^ 
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The Feynman amplitude for the vector boson self-energy graph with two vertices is more com- 
plicated: 

J (27r)3 [/c2 - 1 + ie] [{p - k)^ - + ie] 

Again, there is a factor of three due to the three possible polarizations of the vector bosons on 
the external legs but the important fact is that the tensorial amplitude can be written as the 
sum of three integrals: 

d?k mk^'k" f d^k eiipf'k" + k^p") 



+ 



(27r)3 [k^ - 1 + ie][{p - ky - k'^ + ie] J (27r)3 [{p - k)'^ - + ie][k^ - I + ie] 
d^k Sipf'p" 



(27r)3 [{p - - k2 + ie] [fc2 - 1 + ie] 
Ii'^{p) is in turn the sum of two integrals: 



+ 



+ 



d^k 




Uik/'k" 




(2^)3 

d^k 


[A;2 


- 1 + ie][{p - k)"^ - 1 + ie] 

Uiki^k^iK^ - 1) 




(2^)3 

d^k 


[A;2 


-1 + ie] [{p - A;)2 - 1 + ie] [{p - 
Aik'^g'"' 


ky — K? + ie] 


(2^)3 

d^k 


[A;2 


- 1 + ie][{p - k)"^ - 1 + ie] 

Aik'^g'"'{K^ - 1) 




(27r)3 


[A;2 


-1 + ie] [{p -ky -1 + ie] [{p - 


/c)2 — k2 _|_ ie] 



where the identity k^^k'^ = ^k'^g'^'^ , valid under the integral symbol, has been used. The second 
integral is convergent but the first integral can be written as: 

d^k Aik'^g^"' f d^k Aigi"" 



(27r)3 [k"^ - I + ie][{p - ky - I + ie] J (27r)3 [k'^-l + ie] 

d^k Aig'"'{p^ - 2pk - 1 + ie) 



(27r)3 [A;2 - 1 + ie] [{p - ky - 1 + ie] 
Therefore, 

= Agf" ■ 1(1) + finite part 

because the second integral above can be shown to be finite using the Feyman parametrization 
as in previous subsections. 

To study the li^'^ {p) integral we use Feynman parametrization: 



d^k -Giip'^k" + k/'p" 



(27r)3 [{p -ky-K^ + ie] [k^-l + ie] 
d^k -Qijp^k" + fc^p'^) 
(27r)3 [{k - xpy + 11^ + iey 



1 [ d^k -Giip^^k'^ + k^^pn 2 2 , 2 



Therefore, 

, , f d^q -Giipf'q" + q^p'' + 2xp^'p'') , f d^q -Uixpf'p 

I2 (P) = dx ——^ r o , ,.2 , „-^i2 = dx • 



J (27r)3 [g2 + ^2 + ^£]2 - y (27r)3 [g2 + ^2 + i£]2 

is finite. Here, we have taken into account that integrals with an odd integrand in q^^ are zero. 
^s'^ip) is obviously convergent. 

Adding the three summands, the result used in subsection §. 7.3 is obtained: 2ig^^'^ ■ 7(1) -|- 
finite part. 
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